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1.0  Research  Description 


The  main  objective  of  this  research  is  to  continue  our  theoretical  and  experimental 
investigation  on  physics  of  self-pumped  and  mutually-pumped  phase  conjugators  and  the 
development  of  novel  applications.  In  addition,  we  also  expand  the  scope  of  our  research 
to  include  the  effect  of  partial  coherence  (finite  coherence  length)  on  nonlinear  wave 
mixing  and  phase  conjugation  in  photorefractive  media.  This  is  an  important  and  practical 
subject  which  has  not  been  investigated  by  previous  researchers.  Specifically  we  are 
interested  in  the  effect  of  finite  mutual  coherence  on  the  formation  of  index  gratings, 
especially  2k-gratings,  in  two-wave  mixing  and  self-pumped  phase  conjugation.  The 
same  formulation  and  results  can  be  applied  to  stimulated  Brillouin  scattering  (SBS)  in 
Kerr  media.  The  original  proposal  requested  a  performance  period  of  36  months.  Due  to 
a  budget  cut  at  the  AFOSR,  our  proposal  was  reduced  to  a  8-month  program  with  a 
significantly  lower  funding  level.  During  the  research  period  only  8  months,  we  have 
carried  out  intensive  theoretical  investigation  in  several  important  areas  in  nonlinear 
optical  wave  mixing  and  phase  conjugation,  especially  the  problem  of  beam  coherence. 
Specifically,  we  have  investigated  the  following  problems: 

•  Contradirectional  two-wave  mixing  with  partially  coherent  waves 

•  Two-wave  mixing  with  partially  coherent  waves  in  high-speed  photorefractive 

media 

•  Effect  of  partial  coherence  on  phase  conjugation 

1 . 1  Scientific  Problem 

The  unique  distortion  correction  property  of  optical  phase  conjugation  makes  it 
extremely  attractive  and  important  in  many  advanced  optical  systems.  These  include  free 
space  optical  communication,  atmospheric  correction  and  beam  control  of  high  power 
lasers,  optical  image  processing,  etc.  There  are  many  approaches  which  can  be  employed 
to  achieve  nonlinear  optical  phase  conjugation.  These  include  optical  four-wave  mixing 
(FWM),  stimulated  Brillouin  scattering  (SBS),  self-pumped  phase  conjugation  (SPPC), 
etc.  In  nonlinear  media  such  as  photorefractive  crystals  and  Kerr  media,  it  is  possible  to 
achieve  self-pumped  phase  conjugation  (SPPC)  and  mutually  pumped  phase  conjugation 
(MPPC)[l,2j.  Such  phase  conjugators  are  ideal  for  optical  communications  through 
scattering  media  such  as  atmosphere  or  sea  water.  MPPC  can  also  be  employed  for  the 
phase  locking  of  an  array  of  lasers. 

Both  SPPC  and  MPPC  in  photorefractive  media  have  been  investigated  both 
theoretically  and  experimentally.  Only  transmission  gratings  were  considered  in  the 
theoretical  investigations  for  reasons  of  mathematical  simplicity.  However,  in  most 
experimental  situations,  both  reflection  gratings  and  2k-gratings  are  present  even  when 
the  coherent  length  of  the  laser  is  limited.  A  complete  theory  requires  the  inclusion  of 
reflection  gratings  and  2k-gratings  in  the  coupled  equations.  During  the  current  research 
program,  we  have  investigated  the  role  of  2k-gratings  and  reflection  gratings  in 
photorefractive  mutually  pumped  phase  conjugators  (MPPC).  We  discovered  that, 
depending  on  the  crystal  orientation  and  the  sign  of  the  coupling  constant,  the  presence  of 
2k-gratings  and  reflection  gratings  may  enhance  the  phase  conjugation  efficiency  [3,4],  In 
addition,  the  presence  of  these  additional  gratings  may  lead  to  temporal  instability  in 
phase  conjugation  due  to  grating  competition. 

Although  reflection  gratings  and  2k-gratings  can  be  neglected  in  some  situations, 
they  are  significant  in  many  important  applications  including  phase  locking  of  lasers 
using  MPPC.  In  a  preliminary  model,  we  assume  that  the  grating  strength  is  directly 


related  to  the  normalized  mutual  coherence  of  the  input  beams.  Thus,  by  properly 
arranging  the  path  difference  between  the  input  beams,  it  is  possible  to  control  the  mutual 
coherence,  and  thus  the  strength  of  the  reflection  gratings  and  2k-gratings.  For  example,  it 
is  possible  to  use  two  mutually  incoherent  beams  in  MPPC  in  which  both  reflection 
gratings  and  2k-gratings  can  be  neglected.  This  is  not  entirely  correct  in  SPPC.  As  a 
result  of  a  preliminary  study,  we  found  that  the  mutual  coherence  of  the  two  beams 
actually  depends  on  the  coupling  strength.  In  other  words,  the  mutual  coherence  can  be 
enhanced  or  reduced  by  the  coupling.  This  can  lead  to  an  effective  interaction  length 
(length  of  reflection  gratings  and/or  2k-gratings)  which  can  be  very  different  from  the 
coherent  length. 

1.2  Scientific  and  Technical  Approach 

Nonlinear  optical  phase  conjugation  offers  unique  capabilities  and  possibilities  for 
many  new  applications.  These  include  optical  communication  through  scattering  media, 
wavefront  matched  heterodyne  detection,  wavefront  correction  and  beam  control  of  high 
power  laser  beams  (e.g.,  ABL),  optical  image  processing,  etc.  Although  there  have  been 
numerous  theoretical  and  experimental  works  in  this  area,  technical  issues  such  as 
reflection  gratings,  2k-gratings,  dynamic  range  and  stability  of  the  conjugators  remain 
unresolved.  In  the  following,  we  briefly  describe  the  effect  of  2k-grating  in  self-pumped 
phase  conjugators.  We  then  describe  a  theoretical  model  of  the  role  of  2k-gratings  in 
SPPC  and  a  general  formulation  on  the  effect  of  mutual  coherence  on  the  strength  of  2k- 
gratings  in  photorefractive  media. 

In  self-pumped  phase  conjugators  (SPPCs)  the  counterpropagating  beams  needed 
in  a  four-wave  mixing  process  are  generated  by  the  incident  beam  itself.  By  a  process  that 
is  not  entirely  understood,  the  self-generated  pump  beams  are  often  formed  in  a 
configuration  tending  to  optimize  the  phase  conjugated  beam.  SPPCs  have  been 
demonstrated  in  a  variety  of  configurations  [2].  In  SPPC  involving  no  external  mirrors, 
there  are  at  least  three  different  models  available  to  explain  the  origin  of  the  phase 
conjugation  in  a  single  photorefractive  crystal.  These  are  four-wave  mixing  with  total 
internal  reflection  (FWM-TIR)  [5,6],  stimulated  photorefractive  backscattering  (SPB)  [7- 
9],  four  wave  mixing  and  stimulated  photorefractive  backscattering  (FWM-SPB)  [10-14], 
In  the  last  model,  the  SPPC  formation  relies  on  both  four  wave  mixing  (FWM)  and 
stimulated  photorefractive  backscattering  (SPB).  The  phase  conjugate  beam  is  generated 
by  a  four-wave  mixing  process  involving  the  incident  beam,  forward-propagating  beam 
(the  fanning  beam)  and  its  backward-stimulated  scattering  beam.  In  a  manner  very  similar 
to  stimulated  Brillouin  scattering  (SBS),  the  backward-stimulated  scattering  beam  is 
generated  by  the  SPB  process  involving  2k-gratings.  Hence,  instead  of  a  closed  loop 
inside  the  crystal,  filaments  representing  counter-propagating  beams  are  often  observed. 

In  a  typical  SPPC  experimental  configuration,  an  extraordinarily-polarized  input 
beam  is  incident  on  the  a-face  of  a  regular-cut  BaTi(>3  photorefractive  crystal.  The  SPPC 
formation  mechanism  is  usually  determined  by  the  fanning  pattern  and  the  boundary 
conditions[12].  Generally  speaking,  SPPC  is  often  generated  via  the  SPB  or  FWM-SPB 
mechanisms  in  some  doped  crystals,  but  via  the  FWM-TIR  mechanism  in  undoped 
crystals.  This  is  mainly  due  to  a  significantly  stronger  fanning  and  larger  coupling 
coefficients  for  2k-gratings  in  doped  crystals.  Changing  the  dopant  concentration  or  the 
operating  wavelength  can  cause  a  reconfiguration  of  the  initial  fanning  pattern  which 
leads  to  a  change  of  the  physical  mechanism  and  the  reflectivity  of  SPPC  [12-15]. 
Depending  on  the  boundary  conditions,  the  fanned  beam  may  or  may  not  reach  the  comer 
of  the  crystal  cube.  SPPC  via  FWM-TIR  often  involves  the  participation  of  the  crystal 
comer  to  provide  the  total  reflection. 


In  the  previous  studies,  we  assumed  that  the  two  beams  responsible  for  the  2k- 
gratings  are  mutually  coherent.  This  is  a  good  assumption  provided  that  the  coherence 
length  of  the  laser  is  much  larger  than  the  crystal  dimension.  In  the  FWM-SPB  mode  of 
SPPC,  the  two  beams  are  almost  always  coherent.  This  is  not  entirely  true  in  some 
MPPCs  and  SPPCs.  In  the  event  when  the  coherent  length  of  the  laser  is  similar  to  (or 
even  smaller  than)  the  crystal  dimension,  the  strength  of  the  2k-gratings  will  be  reduced 
due  to  the  reduction  in  the  modulation  depth  (fringe  contrast)  of  the  interference  pattern. 
Assuming  no  beam  coupling,  it  is  known  that  the  modulation  depth  (interference  fringe 
contrast)  is  proportional  to  the  normalized  mutual  coherence  of  the  two  beams.  In  a 
photorefractive  medium,  the  photo-induced  index  grating  is  proportional  to  the 
modulation  depth  (fringe  contrast).  Thus,  it  is  nature  to  assume  that  the  beam  coupling 
constant  is  directly  proportional  to  the  mutual  coherence  of  the  two  beams.  In  other 
words. 


y~l  ~  7^20^12  (2)  (1) 

where  y20  is  the  coupling  constant  for  monochromatic  beams,  T12(z)  is  the  mutual 
coherence  of  the  two  beams  which  is  a  function  of  position  z.  In  a  contradirectional  two- 
wave  mixing,  T12(z)  is  significant  only  over  a  range  of  the  coherence  length.  Thus,  it  is 
often  assumed  that  the  interaction  length  is  limited  to  the  coherence  length.  This 
approximation  is  only  valid  when  the  mutual  coherence  remains  unchanged  during  the 
interaction.  In  the  following  section,  we  will  show  that  the  mutual  coherence  of  the  two 
interacting  beams  can  be  significantly  changed  as  a  result  of  the  beam  coupling, 
especially  when  the  coupling  constant  is  large. 

Two-Beam  Coupling  without  pump  depletion 

Two-wave  mixing  in  photorefractive  crystals  has  been  investigated  extensively  by 
previous  workers  [1].  Most  of  the  theoretical  works  in  this  area  are  on  the  interaction  of 
two  monochromatic  beams.  Two  wave  mixing  with  partially  coherent  beams  has  been 
studied  recently  for  the  case  of  transmission  grating  interaction  only  [19],  for  reasons  of 
mathematical  simplicity.  In  the  case  of  transmission  grating  interaction,  the  optical  path 
difference  between  the  two  interacting  waves  remains  approximately  the  same  as  the  two 
waves  propagate  codirectionally  through  the  photorefractive  medium,  especially  when  the 
incident  angles  of  the  two  waves  are  close  to  each  other.  In  the  case  of  reflection  grating 
interaction,  the  optical  path  difference  between  the  two  interacting  waves  varies 
significantly  as  the  two  waves  propagate  contradirectionally  through  the  photorefractive 
medium.  This  is  the  reason  why  the  mutual  coherence  T12(z)  depends  on  position.  Thus, 
for  the  case  of  transmission  grating  interaction,  only  one  free  variable  for  the  position  is 
needed  to  describe  in  a  self-consistent  way  the  second  order  statistical  properties  of  the 
two  optical  waves,  i.e.  their  intensities  and  mutual  coherence,  while  at  least  two  free 
variables,  one  for  the  position  and  one  for  the  optical  path  difference,  will  be  needed  for 
the  case  of  reflection  grating  interaction.  Another  difficulty  in  studying  the  reflection 
grating  interaction  of  partially  coherent  waves  is  to  find  a  way  to  incorporate  the 
complete  boundary  conditions  into  the  theoretical  formulation  as  a  result  of  the  two-point 
boundary  value  problem.  In  such  a  problem,  a  complete  set  of  boundary  conditions 
includes  both  the  second  order  self  statistical  properties  (e.g.  self-coherence)  of  each 
wave  at  its  entrance  boundary  and  the  second  order  mutual  statistical  properties  (e.g. 
mutual  coherence)  of  the  two  waves  at  their  respective  entrance  boundaries.  In  this 
section,  we  present  a  preliminary  theoretical  investigation  for  the  case  of  reflection 
grating  interaction.  In  the  case  of  contradirectional  two  wave  mixing  with  purely 


diffusion-driven  photorefractive  medium,  the  coupled-wave  equations  for  the  slowly 
varying  amplitudes  E,(t,z)  and  E2(t,z)  can  be  written  as  [19] 


3E 
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dz 

v  9t 

2I,  +  I2 

3E2 

1  3E2 
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where  y  is  the  intensity  coupling  constant,  a  is  the  intensity  absorption  coefficient,  v  is 
the  group  velocity,  E^xpHcot  +  ikz)  and  E2exp(-icot  -  ikz)  are  the  coupled  quasi- 
monochromatic  waves  which  interact  through  a  dynamic  photorefractive  grating 
8n(z,t)  oe  Qexp[2ikz]  +  c.c..  The  dynamics  of  the  photorefractive  index  grating  can  be 
written  as 


t^  +  Q  =  E,E2*  (4) 

3t 

where  t  is  the  relaxation  time  constant  determined  by  the  average  total  intensity, 
T~  ((£;)2  +(E2)2) 1 .  Reflecting  the  nature  of  partial  coherence,  the  complex  amplitudes 
of  the  two  waves  are  considered  random  variables,  especially  in  the  regime  when  the 
coherence  time  8o)_1  is  substantially  less  than  the  relaxation  time  of  the  material 
Sore » 1  [20],  When  the  optical  path  difference  of  the  two  beams  is  smaller  than  the 
coherence  length  of  the  laser,  a  dynamic  photorefractive  grating  is  recorded  in  the 
medium.  Its  position  and  profile  are  nearly  stationary.  Thus,  as  an  approximation,  we  can 
replace  the  dynamic  grating  amplitude  Q  in  Eqs.(2)  and  (3)  with  its  ensemble  average  as 

Q-  (E,E,-). 

Since  the  complex  amplitudes  E^z)  and  E2(t,z)  are  stationary  random 
processes,  we  can  define  some  of  their  ensemble  averages  as 
ri2  (z,  At)  s  (E,  (z,  t,  )E2*  (z,  t2 )),  f,  i  (z,  At)  s  (e,  (z,  t,  )Ej*  (z,  t2 )) ,  and 

r22(z,At)H(E2(z,t1)E2*(z,t2))  where  At  =  t,  - 12  and  (  )  means  ensemble  average.  We 
refer  to  ru  and  T22  as  the  self  coherence  of  the  pump  beam  E,(t,z)  and  the  signal  beam 
E2(t,z)  respectively,  and  T12  as  the  mutual  coherence  between  the  signal  beam  and  the 
pump  beam.  With  these  notations,  we  can  immediately  write  down  Q  =  ri2(z,0), 
lj(z)  =  rn(z,0)  and  I2(z)  =  r22(z,0).  Eqs.(2)  and  (3)  can  therefore  be  reduced  to  a 
system  for  these  averaged  values  as 


£^^=xL^rn-(z,-At)+|if^)r1!(z,At)-arll(z,At)  (6) 
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(z,  At)  +  af22  (z,  At)  (7) 


Note  that  Eqs.(5)-(7)  represent  the  propagation  of  the  mutual  conference  and  self 
coherence  functions  of  the  two  beams.  When  there  is  no  coupling  (i.e.  y  =  0),  Eqs.(5)-(7) 
can  be  solved  analytically.  The  solutions  are  consistent  with  the  classical  theory  on 
optical  wave  interference  in  linear  medium[21].  The  mutual  coherence  function  is 
presented  as  a  propagation  wave  with  the  boundary  conditions  at  z  =  0  plane.  The  self 
coherence  function  decreases  along  the  z  axis  due  to  absorption.  When  the  coupling  is 
present,  both  beams  are  scattered  by  the  grating  into  each  other.  As  a  result,  part  of  pump 
wave  branches  off  in  the  direction  of  the  signal  waves.  Thus,  the  mutual  and  self 
conference  functions  of  the  two  waves  can  not  be  derived  by  the  simple  propagation  of 
the  boundary  conditions.  However,  we  note  that  Eqs.(5)-(7)  are  self-consistent.  They  can 
be  solved  numerically  if  we  have  the  proper  boundary  conditions.  In  the  presence  of 
pump  depletion,  we  can  not  obtain  sufficient  boundary  conditions  for  typical  situations. 
Therefore,  although  Eqs.(5)-(7)  are  general,  they  are  not  useful  for  the  case  of  pump 
depletion.  For  the  case  of  nondepleted  pump,  solutions  with  sufficient  boundary 
conditions  can  be  easily  obtained.  Preliminary  results  indicate  that  the  contradirectional 
two  wave  mixing  of  partial  coherent  beams  inside  a  photorefractive  crystal  with  an 
inertial  nonlinear  response  can  lead  to  a  increase  of  the  mutual  coherence  substantially 
and  the  interaction  length  of  the  two  beams.  A  set  of  the  coupled  equation  has  been 
derived  to  describe  the  propagation  of  the  mutual  coherence  and  the  self  coherence 
functions  of  the  two  beams. 


Two-Beam  Coupling  with  pump  depletion 


When  beam  depletion  is  not  negligible,  the  investigation  requires  the  use  of 
statistical  optical  techniques,  especially  for  the  case  of  reflection  gratings.  In  the  above, 
we  describe  a  theoretical  formulation  of  the  problem  in  the  space  and  time  domain  for  the 
reflection  grating  interaction  in  the  nondepleted  pump  regime.  By  using  the  nondepleted 
pump  approximation,  we  reduced  the  two-point  boundary  value  problem  to  an  initial 
value  problem.  In  this  section,  we  present  a  preliminary  result  of  our  general  formulation 
of  the  problem  in  the  space  and  frequency  domain  based  on  the  standard  statistical  theory 
on  linear  systems.  Specifically,  we  are  interested  in  the  signal  intensity  gain  and  the 
mutual  coherence  in  the  contradirectional  wave  mixing  of  two  partially  coherent  waves. 
Photorefractive  two- wave  mixing  is  a  nonlinear  optical  process.  Because  of  the  mutual 
coherence  of  the  two  waves,  a  dynamic  holographic  grating  is  formed  in  the  medium.  Its 
position  and  index  profile  are  nearly  stationary  under  the  condition  of  a  cw  illumination. 
Both  waves  are  scattered  into  each  other  by  the  presence  of  this  index  grating.  Scattering 
of  partially  coherent  waves  by  a  stationary  grating  can  be  modeled  with  a  statistical 
theory  on  linear  systems  [21].  An  iterative  procedure  can  subsequently  be  devised  to 
obtain  the  final  photorefractive  grating  profile  from  an  initially  arbitrary  grating  profile. 


We  consider  the  interaction  of  two  counter-propagating  waves  with  partial 
coherence  in  a  photorefractive  medium  between  z=0  and  z=L.  Reflecting  the  partial 
coherence,  E,(z,t)  and  E2(z,t)  are  treated  as  stationary  random  variables.  They 


represent  the  random  fluctuation  of  the  amplitudes  of  the  two  waves.  Let  Eu(z,Aco)  and 
E22(z,Aco)  denote  the  self  spectral  density  functions  of  Ej(z,t)  and  E2(z,t)  respectively 

and  E12(z!,z2,Aco)  be  the  cross  spectral  density  function  between  E,(z,,t)  and  E2(z2,t). 
The  spectral  density  functions  and  the  corresponding  coherence  functions  are  Fourier 
transform  pairs,  i.e. 


ru(z,x)  =  jEn(z,AcD)exp(-iAcDx)dAG) 

(8) 

r22(z,x)  =  jE22(z,Aco)exp(-iAG>x)dAco 

(9) 

T12(z,x)  =  J  E12(z,z,  Aw)exp(-iA(Ox)dA(D 

(10) 

where  Aco  =  (0  -  co0  with  CD  being  a  general  frequency  component  of  the  waves.  With 
the  above  notations  and  relations,  the  intensity  of  the  two  waves  can  be  expressed  as 

I,(z)  s  rn(z,0)  =  jEu(z,Aco)dAco  (11) 

I2  (z)  h  r22  (z,  0)  =  J  E^  (z,  Ato)dAcD  (12) 

and  the  mutual  coherence  of  the  two  waves  can  be  expressed  as 

T12(z,0)  =  jE12(z,z,Aco)dAcD  (13) 


Using  these  equations,  the  intensity  of  the  waves  as  well  as  the  mutual  coherence  can  be 
obtained  as  soon  as  the  spectral  density  functions  are  obtained.  In  what  follows,  we  will 
derive  the  spectral  density  functions  by  using  the  statistical  approach. 


Using  the  above  notations  and  definitions,  we  now  begin  our  discussion  on  the 
photorefractive  interaction.  Inside  the  photorefractive  medium,  a  dynamic  index  grating  is 
generated.  It  can  be  written  as 


s  .Y  c 

on  =  -i - 

2  con 


^Y7exp(2ik0z)  +  c.c. 

Io(*) 


(14) 


where  Q(z,t)  is  a  measure  of  the  index  grating,  y  is  the  intensity  coupling  coefficient 
and  I0(z)  =  I,(z)  +  I2(z)  is  the  total  intensity  at  position  z.  For  the  purpose  of  our 
discussion,  we  will  call  E,  as  the  signal  wave  and  E2  as  the  pump  wave.  Thus,  for  a 
photorefractive  grating  with  a  positive  y ,  the  signal  wave  Ej  can  be  amplified.  If  the 


photorefractive  effect  is  based  purely  on  carrier  diffusion  (e.g.  BaTi03),  the  dynamics  of 
the  index  grating  is  described  by  Eq.  (3).  By  virtue  of  photoexcitations,  photorefractive 
processes  are  usually  slow  at  low  intensities.  It  is  reasonable  to  assume  that  the  coherence 

time  5co-1  of  the  two  partially  coherent  waves  is  much  smaller  than  the  relaxation  time  r 
of  the  photorefractive  medium,  i.e,  8(0X»  1.  Since  Ej(z,t)  and  E2(z,t)  are  stationary 
random  variables,  we  can  make  the  following  approximation,  according  to  Eq.(3), 


Q(z,t)  =  (Q(z,t))  =  T12(z,0) 


(15) 


In  other  words,  two  partially  coherent  waves  with  their  complex  amplitudes  fluctuating 
randomly  with  time  can  actually  write  a  stationary  grating  in  a  photorefractive  medium 
under  the  appropriate  conditions.  For  simplicity,  we  will  denote  Q(z,t)  with  Q(z)  from 
now  on.  Note  that  Q(z)  is  also  the  mutual  coherence  of  the  two  waves  at  position  z. 
Given  arbitrary  functions  of  Q(z)  and  I0(z),  Eq.(14)  yields  an  index  grating.  The 
propagation  of  a  monochromatic  wave  through  such  an  index  grating  can  be  described  by 
the  coupled  wave  equations.  Let  the  general  solutions  be  written  Hij(z,o))  (i,j=l,2).  We 

note  that  each  input  wave  generates  a  transmitted  and  a  reflected  wave.  The  solutions  can 
then  be  employed  to  modify  the  index  grating  as  a  result  of  the  photorefractive  effect. 

This  leads  to  a  new  index  profile  Q(z) .  The  iterations  continue  until  the  index  profile 
converges. 

In  each  of  the  iterations,  an  arbitrary  stationary  index  grating  can  be  considered  as 
a  linear  system.  We  consider  a  given  stationary  index  grating  in  the  photorefractive 
medium  as  a  linear  system  with  the  optical  waves  at  the  boundary  planes  z=0  and  z=L  as 
the  input  and  the  optical  waves  at  an  arbitrary  plane  z  as  the  output.  According  to  the 
theory  on  the  statistical  properties  of  linear  systems,  the  second  order  statistical  properties 
of  the  optical  waves  at  the  output  plane  can  be  expressed  in  terms  of  the  second  order 
statistical  properties  of  the  optical  waves  at  the  input  planes  and  the  frequency  response 
of  the  linear  system.  Two  physical  processes  happen  simultaneously  during  two-wave 
mixing  in  a  photorefractive  medium.  First,  the  two  optical  waves  propagate  through  the 
photorefractive  medium  while  being  scattered  into  each  other  by  the  index  grating. 
Second,  the  scattered  waves  modify  the  index  grating  through  the  photorefractive  effect 
until  a  steady  state  is  reached.  We  have  provided  above  a  mathematical  model  that 
describes  these  two  physical  processes  separately.  A  steady  state  of  the  two-wave  mixing 
in  the  photorefractive  medium  is  reached  when  the  two  optical  waves  scattered  by  the 
photorefractive  grating  can  exactly  sustain  the  same  photorefractive  grating.  A  steady 
state  solution  of  the  two-wave  mixing  in  a  photorefractive  medium  can  thus  be  obtained 
by  using  the  mathematical  model  described  above  through  an  iterative  procedure.  The 
research  results  in  this  program  are  based  on  the  statistical  approach  described  above. 

1.3  Publications: 

1.3.1  Papers  Published: 


X.  Yi,  C.  Yang,  P.  Yeh,  S.  Lin,  K.  Y.  Hsu,  "General  solution  of 
contradirectional  two-wave  mixing  with  partially  coherent  waves  in 


photorefractive  crystals, "  J.  Opt.  Soc.  Am.,  B,  Vol.  14,  No.  6, 1396-1406 
(1997). 

X.  Yi  and  P.  Yeh,  "Two-wave  mixing  with  partially  coherent  waves  in  high¬ 
speed  photorefractive  media, "  J.  Opt.  Soc.  Am.,  B,  Vol.  14,  No.  11,  2885-2894 
(1997). 

X.  Yi  and  P.  Yeh,  "Effect  of  partial  coherence  on  phase  conjuagtion,"  Opt. 
Comm.,  Vol.  147,  Nl-3, 126-130  (1998). 
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X.  Yi  and  P.  Yeh,  "Phase  conjugation  through  turbulent  atmosphere,"  1997 
OS  A  Annual  Meeting  (Long  Beach,  CA,  October  13-17, 1997),  paper  MJJ4. 

R.  Wang  and  P.  Yeh,  "Optical  power  limiting  via  both  two-wave  mixing  and 
beam  fanning  in  BaTi03  crystals,"  1997  OSA  Annual  Meeting  (Long  Beach, 
CA,  October  13-17, 1997),  paper  ThFF24. 

X.  Yi,  P.  Yeh  and  S.  Campbell,  "Effect  of  thermal  expansion  on  holographic 
memory,"  1997  OSA  Annual  Meeting  (Long  Beach,  CA,  October  13-17,  1997), 
paper  FK4. 


P.  Yeh,  "Nonlinear  Optical  Wave  Mixing  with  Partially  Coherent  Waves," 
(Invited  Paper)  Waseda  International  Symposium  on  Phase  Conjugation  and 
Wave  Mixing,  June  9-10,  1997  (Waseda  University,  Tokyo,  Japan). 

P.  Yeh,  "Recent  Advances  in  Photorefractive  Nonlinear  Optics,"  (Invited 
Paper)  CLEO/Pacific  Rim'97  (July  14-18,  Tokyo,  Japan,  1997). 


2.0  Progress: 

2. 1  Progress  Summary: 

2.1.1  Contradirectional  two- wave  mixing  with  partially  coherent  waves 

During  the  period  we  investigated  the  contradirectional  two-wave  mixing  with  partially 
coherent  waves  in  photorefractive  media.  We  first  considered  the  case  when  the  material 
response  is  slow.  This  is  typical  of  the  photorefractive  crystals  with  typical  laser  powers. 
By  use  of  the  statistical  method  described  above,  a  general  solution  of  the  problem  in  the 
space  and  frequency  domain  is  obtained.  We  obtained  results  on  beam  intensity  and 
mutual  coherence.  The  results  on  the  enhancement  of  beam  coherence  are  compared  with 
previous  theoretical  results  on  simpler  cases  and  with  experimental  measurements. 
Excellent  agreements  are  achieved.  The  results  also  indicate  that  the  effective  interaction 
length  can  be  significantly  longer  than  the  coherence  length  of  the  waves. 

2.1.2  Two- wave  mixing  in  high  speed  media 

In  the  case  when  the  laser  power  is  very  high,  the  material  response  is  proportionally 
fast.  We  investigated  the  coupling  of  two  partially  coherent  waves  in  photorefractive 


media  with  a  fast  temporal  response.  When  the  coherence  time  of  the  waves  is  close  to 
the  speed  of  the  photorefractive  media,  the  photo-induced  index  grating  will  follow  the 
amplitude  fluctuation  of  the  optical  waves.  The  stochastic  coupled  wave  equations  are 
solved  in  the  non-depleted  pump  region.  The  effects  on  the  evolution  of  the  signal  wave 
are  obtained. 

2.1.3  Effect  of  partial  coherence  on  phase  conjugation 

During  the  period,  we  also  investigated  the  effect  of  partial  temporal  coherence  of  the 
incident  wave  on  the  phase  conjugation  process  based  on  2k-grating.  In  many 
photorefractive  phase  conjugators,  2k-gratings  are  generated  by  the  incident  beam  and  the 
phase  conjugated  beam.  The  phase  conjugated  beam  will  be  amplified  by  the  incident 
beam  via  contra-directional  two-wave  mixing.  The  strength  of  the  2k-grating  and 
therefore  the  phase  conjugate  reflectivity  will  be  dependent  on  the  coupling  between 
them.  We  obtained  both  theoretical  and  experimental  results  of  the  effect  of  beam 
coherence  on  the  phase  conjugate  reflectivity. 

2.2  Progress  details: 

Since  most  of  the  details  of  research  results  have  already  been  published  in 
technical  journals.  Detail  discussion  on  the  research  results  can  also  be  found  in  Section  4 
which  contains  reprints  of  the  papers  published. 
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We  investigate  the  coupling  of  two  partially  coherent  waves  in  high-speed  photorefractive  media.  When  the 
coherence  time  of  the  optical  waves  is  close  to  the  speed  of  the  photorefractive  media,  the  photoinduced  index 
grating  will  follow  the  amplitude  fluctuation  of  the  optical  waves.  The  stochastic  coupled- wave  equations  are 
solved  in  the  nondepleted  pump  region.  The  effects  on  the  evolution  of  the  signal  wave  are  presented  and 
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1.  INTRODUCTION 

Nonlinear  propagation  and  interaction  of  partially  coher¬ 
ent  optical  waves  in  photorefractive  materials  have  re¬ 
ceived  a  lot  of  attention  in  recent  years,  because  of  appli¬ 
cations  including  image  amplification,  optical 
interconnection,  optical  phase  conjugation,  and  photore¬ 
fractive  filters.1"6  In  most  previous  theoretical  analysis, 
photorefractive  materials  have  been  assumed  to  have  a 
long  response  time  r  so  that  their  photorefractive  nonlin¬ 
earity  does  not  follow  the  instantaneous  fluctuation  of  the 
partially  coherent  optical  waves.7-15  However,  the  re¬ 
sponse  time  of  the  photorefractive  materials  depends  both 
on  the  materials  themselves  and  on  the  intensities  of  the 
interacting  waves.16  Furthermore,  whether  the  instan¬ 
taneous  photorefractive  nonlinear  response  can  be  ne¬ 
glected  depends  on  the  relative  photorefractive  response 
time  with  respect  to  the  bandwidth  So)  of  the  interacting 
waves  rather  than  on  the  photorefractive  response  time 
itself.  The  assumption  of  neglecting  the  instantaneous 
photorefractive  response  is  valid  mainly  for  low  optical  in¬ 
tensities  and  wide  optical  bandwidths  so  that  the  condi¬ 
tion  Sa)T  >  1  is  satisfied.  When  the  optical  intensity  is 
high  and  the  optical  bandwidth  is  narrow  so  that  Swr  is 
of  the  order  of  unity,  the  photorefractive  nonlinearity  will 
respond  to  the  instantaneous  fluctuation  of  the  partially 
coherent  optical  waves.  The  nonlinear  optical  properties 
in  the  high-speed  case  are  significantly  different  from 
those  in  the  low-speed  case. 

In  this  paper  we  investigate  codirectional  two-wave 
mixing  with  partially  coherent  waves  in  high-speed  pho¬ 
torefractive  media.  When  the  spectral  bandwidth  of  the 
optical  waves  is  close  to  the  speed  of  the  photorefractive 
media,  the  photoinduced  index  grating  will  fluctuate  with 
the  partially  coherent  optical  waves.  This  case  is  much 
more  complicated  than  the  corresponding  case  in  low- 
speed  photorefractive  media.  To  simplify  the  problem, 
we  limit  our  analysis  to  the  nondepleted  pump  region  and 
assume  that  the  pump  wave  has  only  phase  noise.  The 
effects  on  the  evolution  of  the  signal  wave  are  determined. 
The  results  show  that  beam  coupling  can  still  be  achieved 
for  incoherent  beams  provided  that  the  material  reacts 
fast  enough.  It  is  important  to  distinguish  wave  mixing 
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due  to  shared  noise  gratings9"11  from  the  direct  coupling 
of  incoherent  beams  described  here. 

This  paper  is  divided  into  three  sections.  In  Section  2, 
a  set  of  stochastic  coupled-wave  equations  is  introduced 
along  with  the  necessary  assumptions  that  ensure  their 
validity  in  the  high-speed  case.17  In  Section  3,  the  sto¬ 
chastic  coupled-wave  equations  are  reduced  to  nonsto¬ 
chastic  differential  equations,  from  which  the  signal  in¬ 
tensity  and  the  mutual  coherence  of  the  two  waves  are 
obtained.  In  Section  4,  the  stochastic  coupled-wave 
equations  are  rewritten  into  a  series  form,  from  which  an 
asymptotic  solution  of  the  signal  wave  bandwidth  is  ob¬ 
tained. 


2.  THEORETICAL  MODEL 

In  this  section  we  present  a  theoretical  model  of  wave 
mixing  in  the  medium.  Several  major  assumptions  are 
made  and  discussed  as  we  derive  the  theoretical  model. 

The  system  configuration  is  shown  in  Fig.  1.  The  sig¬ 
nal  wave  and  the  pump  wave  enter  the  photorefractive 
medium  at  the  incidence  plane  2  =  0.  The  bisector  of  the 
two  beams  is  perpendicular  to  the  surface  normal.  A 
transmission  grating  with  grating  wave  vector  along  the  x 
direction  is  generated  inside  the  photorefractive  medium. 
The  two  optical  waves  are  coupled  through  the  index  grat¬ 
ing. 

We  assume  that  the  two  optical  waves  are  plane  waves 
of  infinite  extent  in  the  transverse  dimension.  This  is 
the  first  major  assumption  we  make.  Under  this  as¬ 
sumption  the  evolution  of  the  optical-wave  amplitudes 
can  be  confined  to  only  one  dimension.  Thus,  the  optical 
waves  inside  the  photorefractive  medium  can  be  written 
as 

E(z,  t)  =  Ex{z ,  t)exp(-ia)0t  +  -  ikopc) 

+  E2(z ,  t)\exp(~ia)0t  +  ikozz  +  iko**), 

(1) 

where  E(z,  t)  is  the  total  electric  field  of  the  two  waves, 
Ei  is  the  signal  wave  complex  amplitude,  E2  is  the  pump- 
wave  complex  amplitude,  o>0  is  the  center  frequency  of  the 
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Fig.  1.  System  configuration:  codirectional  two-wave  mixing 
with  partially  coherent  waves. 

two  partially  coherent  waves,  ko  =  tkjjq/c  is  the  center- 
wave  vector,  n  is  the  refractive  index  of  the  photorefrac- 
tive  medium,  and  c  is  the  speed  of  light  in  vacuum. 
Reflecting  the  partial  coherence,  Ex(zt  t)  andE2{zf  t)  are 
stationary  random  variables.  They  represent  the  ran¬ 
dom  fluctuation  of  the  amplitudes  of  the  two  waves. 

The  instantaneous  intensities  of  the  two  waves  are 
Ex(z,  t)Ex*(z ,  t)  and  E2(z,  t)E2*(z9  t ).  In  general, 
they  fluctuate  randomly  with  time,  e.g.,  as  in  thermal 
light  and  multimode  laser  light.  The  intensities  of  the 
two  waves  are  defined  as  the  ensemble  average  of  their 
respective  instantaneous  intensities: 


tal  intensity  of  the  two  waves  need  to  be  taken  into  ac¬ 
count  to  allow  us  to  determine  the  photorefractive  grat¬ 
ing.  This  is  a  difficult  task  in  general.  To  overcome  this 
difficulty,  we  introduce  two  more  major  assumptions. 
One  is  that  the  pump  wave  is  much  stronger  than  the  sig¬ 
nal  wave  and  is  undepleted  in  the  course  of  the  interac¬ 
tion.  The  other  is  that  the  pump  wave  has  only  phase 
noise  at  the  incidence  plane.  Under  these  two  assump¬ 
tions,  the  pump  wave  will  have  only  phase  noise  through¬ 
out  the  photorefractive  medium.  The  instantaneous  to¬ 
tal  intensity  of  the  two  waves  is  approximately  equal  to 
the  instantaneous  intensity  of  the  pump  wave  and  is  a 
constant  over  time  throughout  the  photorefractive  me¬ 
dium.  We  now  need  only  to  consider  the  instantaneous 
mutual  coherence  of  the  two  waves  to  determine  the  pho¬ 
torefractive  grating.  The  photorefractive  grating  can  be 
written  as 

y  £  Q(z  i) 

'  Sn(z,  t)  =  -i  - - exp(-2fk0xx)  +  cc,  (5) 

2  wo  I2(z) 

where  Q(z,  t)  is  the  photorefractive  grating  amplitude 
factor.  In  a  photorefractive  medium  of  purely  diffusive 
type,  the  dynamics  of  the  photorefractive  grating  can  be 
written  as 


Ii(z)  -  {E\{z ,  t)EY*(z%  t))9  (2) 

I2(z)  -  (E2(z,  t)E2*(zy  t)).  (3) 

The  intensities  of  the  two  waves  as  defined  by  Eqs.  (2) 
and  (3)  are  always  constant  over  time.  There  is  a  special 
case.  When  a  partially  coherent  optical  wave  has  only 
phase  noise,  its  instantaneous  intensity  is  a  constant  over 
time  and  is  equal  to  its  ensemble  averaged  intensity,  e.g., 
for  single-mode  laser  light.  The  instantaneous  mutual 
coherence  of  the  two  waves  is  Ex(z}  t)E2*(z,  t),  which 
also  fluctuates  randomly  with  time.  The  mutual  coher¬ 
ence  of  the  two  waves  is  defined  as  the  ensemble  average 
of  the  instantaneous  mutual  coherence  of  the  two  waves. 

We  now  consider  the  photorefractive  grating  dynamics. 
In  the  case  of  two-wave  mixing  with  partially  coherent 
waves  in  low-speed  photorefractive  media,  the  photore¬ 
fractive  grating  can  be  written  as 

y  c  ( Ex(z ,  t)E2*(z,  t)) 

Sn{z,  t)  -  i2o)o  /,(*)  +  /,(*) 

x  exp(-2ik0xx)  +  cc,  (4) 

where  y  is  the  intensity  coupling  coefficient.  The  ampli¬ 
tude  of  the  photorefractive  grating  is  proportional  to  the 
ensemble  averaged  mutual  coherence  of  the  two  waves 
and  is  inversely  proportional  to  the  ensemble  averaged  to¬ 
tal  intensity  of  the  two  waves.  In  low-speed  photorefrac¬ 
tive  media  the  photorefractive  nonlinearity  does  not  re¬ 
spond  to  the  instantaneous  fluctuations  of  the  partially 
coherent  optical  waves.  Equation  (4)  is  a  valid  approxi¬ 
mation  for  most  kinds  of  partially  coherent  wave,  e.g., 
thermal  light,  multimode  laser  light,  and  single-mode  la¬ 
ser  light.  In  high-speed  photorefractive  media,  the  pho¬ 
torefractive  nonlinearity  responds  to  the  instantaneous 
fluctuation  of  the  partially  coherent  optical  waves.  Both 
the  fluctuation  of  the  instantaneous  mutual  coherence  of 
the  two  waves  and  the  fluctuation  of  the  instantaneous  to¬ 


t  dQ(Z’  -  +  Q(z,  t)  =  Ei{z,  t)E2*(z,  t),  (6) 

dt 

where  r  is  the  photorefractive  grating  relaxation  time 
constant.  Equation  (6)  can  be  written  in  an  integral 
form, 

Q(z,  t)  =  ^  J  Ey(z,  t')E2*(z,  Oexpj  — — jdf\ 

(7) 

Equation  (7)  is  mathematically  equivalent  to  Eq.  (6)  and 
is  sometimes  more  convenient  to  use  than  Eq.  (6)  in  our 
analysis. 

With  the  photorefractive  grating  given  by  Eq.  (5)  and 
the  undepleted  pump  assumption,  we  can  write  the 
coupled  wave  equations  as 


dEx(zy  t)  1  dEx(zf  t) 
dz  v  dt 


y  Q(z,  t) 
2  l2(z) 


E2(Zy 


t)> 


(8) 


dE2{zy  t)  +  1  dE2(zy  t)  =  o  (9) 

dz  V  dt 

where  v  is  the  group  velocity  of  the  two  waves  in  the  pho¬ 
torefractive  medium.  In  arriving  at  Eqs.  (8)  and  (9),  we 
have  also  assumed  that  there  is  no  absorption  in  the  pho¬ 
torefractive  medium.  This  is  the  fourth  major  assump¬ 
tion  needed  to  simplify  the  analysis.  A  linear  absorption 
can  later  be  incorporated  into  the  coupled  wave  equa¬ 
tions. 

In  the  above  derivation  we  have  made  no  assumption 
about  the  statistical  properties  of  the  signal  wave.  The 
signal  wave  is  allowed  to  have  both  intensity  noise  and 
phase  noise  at  the  incidence  plane.  In  addition,  the  sta¬ 
tistical  properties  of  the  signal  wave  are  also  allowed  to 
change  in  the  course  of  interaction.  The  assumption  that 
the  pump  wave  at  the  incidence  plane  has  only  phase 
noise  can  be  justified  with  typical  experimental  condi- 
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tions.  According  to  Ref.  16,  the  photorefractive  grating 
relaxation  time  constant  decreases  as  the  total  intensity 
of  the  two  interacting  waves  increases.  When  the  total 
intensity  is  low,  e.g.,  1  W/cm2,  the  grating  relaxation  time 
constant  falls  into  the  range  from  1  ms  to  1  s.  When  the 
total  intensity  is  high,  the  grating  relaxation  time  con¬ 
stant  can  be  of  the  order  of  1  /ns  or  lower.  The  bandwidth 
of  an  optical  wave  whose  coherence  time  is  close  to  the 
grating  relaxation  time  at  the  short  end  will  be  of  the  or¬ 
der  of  1  MHz.  This  is  a  narrow  bandwidth  even  for  a 
single-mode  laser  source.  The  assumption  that  the  pump 
wave  at  the  incidence  plane  has  only  phase  noise  agrees 
well  with  experiments  in  which  the  pump  wave  is  ob¬ 
tained  from  a  single-mode  laser.  As  we  shall  see,  this  as¬ 
sumption  and  the  undepleted  pump  assumption  lead  to 
significant  simplifications  for  solving  the  coupled-wave 
equations. 

To  solve  the  coupled-wave  equations,  we  need  to  know 
certain  information  about  the  statistical  properties  of  the 
signal  wave  and  the  pump  wave  at  the  incidence  plane  as 
the  boundary  condition.  Without  loss  of  generality,  we 
assume  that  the  signal  wave  and  the  pump  wave  are  ob¬ 
tained  from  the  same  source  laser  wave.  We  also  assume 
a  single-mode  laser  source  and  that  the  source  laser  wave 
has  only  phase  noise.  The  optical  wave  amplitudes 
E^z,  t)  and  E2(z7  t)  at  the  incidence  plane  can  thus  be 
written 

E^O.t)  =  VA( 0)  exp[-i0(f  +  fd)],  (10) 

E2(0,t)  =  sfh(0)  exp [-*«*)].  (11) 

where  td  is  the  time  delay  between  the  signal  wave  and 
the  pump  wave  at  the  incidence  plane  and  0{t)  is  a  ran¬ 
dom  variable  representing  the  random  phase  fluctuation 
of  the  optical-wave  amplitudes  at  the  incidence  plane. 
The  random  phase  0(t)  is  assumed  to  have  stationary  in¬ 
crements  [i.e.,  0(t  +  St)  -  0(t)  is  statistically  indepen¬ 
dent  of  the  time  t].  Furthermore,  the  increment  over  an 
arbitrarily  given  time  period  Stf  i.e.,  &0(St)  =  0(t  +  St) 
-  0{t)y  is  assumed  to  have  a  Gaussian  distribution.  The 
random  phase  shifts  over  nonoverlapping  time  intervals 
are  assumed  to  be  statistically  independent.  Under 
these  assumptions,  the  source  laser  wave  has  a  Lorentz- 
ian  line  shape,  i.e.,  the  typical  line  shape  of  a  single-mode 
semiconductor  laser.  The  variance  of  the  random  phase 
shift  over  an  arbitrarily  given  time  period  St  is  related  to 
the  bandwidth  of  the  source  laser  wave  Sw  as 

<A  0(St))2  =  So)St.  (12) 

The  above  assumptions  for  the  source  laser  wave  are  typi¬ 
cal  for  a  single-mode  laser  source. 

The  theoretical  model  is  complete  with  the  coupled 
wave  equations  (7)-(9)  and  the  assumptions  for  the  sta¬ 
tistical  properties  of  the  signal  wave  and  the  pump  wave 
at  the  incidence  plane  as  the  boundary  condition. 

3.  SIGNAL  INTENSITY  AND  MUTUAL 
COHERENCE 

In  this  section  we  solve  for  the  signal  intensity  and  the 
mutual  coherence  of  the  two  waves.  The  nonlinear  sto¬ 
chastic  partial  differential  equations  (7)-(9)  are  reduced 


to  nonstochastic  differential  equations,  the  solutions  of 
which  include  the  signal  intensity  and  the  mutual  coher¬ 
ence  of  the  two  waves.  This  section  is  further  divided 
into  three  subsections. 


A.  Solutions 

In  this  subsection  the  coupled-wave  equations  (7)-(9)  are 
solved,  and  the  expressions  for  the  signal  intensity  and 
the  mutual  coherence  of  the  two  waves  are  obtained  along 
with  the  necessary  boundary  conditions. 

We  first  look  at  the  pump  wave.  At  the  incidence 
plane,  the  pump  wave  has  only  phase  noise.  With  Eq. 
(9),  the  amplitude  of  the  pump  wave  inside  the  photore¬ 
fractive  medium  can  be  written  as 

E2(zy  t)  -  E2(Q,t  -  z/v).  (13) 

The  pump  wave  propagates  through  the  photorefractive 
medium  without  changing  its  statistical  properties. 
Therefore  the  pump  wave  has  the  same  intensity  and  has 
only  phase  noise  throughout  the  photorefractive  medium, 
i.e., 

I%(z)  =  /2(0),  (14) 

E2(z,  t)E2*(z ,  t)  =  /2(0).  (15) 

Note  that  the  instantaneous  intensity  E2{z ,  t)E2*(z7  t) 
is  a  constant  equal  to  the  ensemble-averaged  intensity 
when  E2(zt  t)  has  only  phase  fluctuations. 

Second,  we  consider  the  mutual  coherence  of  the  two 
waves.  We  denote  the  instantaneous  mutual  coherence 
of  the  two  waves  as 

F(Z,  t)  m  Ei(zt  t)E2*(z,  t).  (16) 

Using  Eqs.  (7)-(9),  and  (15),  we  obtain 


dF{zy  t)  1  0F(z,  t) 

Bz  V  Bt 


F(z ,  f) 


X  exp|~jd£'.  (17) 

This  equation  governs  the  propagation  of  the  instanta¬ 
neous  mutual  coherence  of  the  two  waves  in  the  photore¬ 
fractive  medium.  Note  that  Eq.  (17)  is  a  linear  equation 
of  the  random  variable  F(z ,  t).  Taking  the  ensemble  av¬ 
erage  of  Eq.  (17),  we  obtain 


(F(z,  <)>  =  <F(0,O>exp(y).  (18) 

( F(z ,  t))  is  the  ensemble-averaged  mutual  coherence  of 
the  two  waves.  According  to  Eq.  (18),  the  mutual  coher¬ 
ence  of  the  two  waves  increases  exponentially  with  the 
position  z . 

We  now  define  a  fourth-order  coherence  function  of  the 
optical-wave  amplitudes  t)  and  E2{zy  t)  as 

R(z,  A t)  m  {Ex(z7  t  +  A t)E2*(zf  t  +  M)Ex*(z9  t) 

x  E2(z ,  t)).  (19) 

According  to  Eqs.  (15)  and  (19),  we  have 


hiz) 


R(Zy  0) 

/2(0) 


(20) 
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Note  that  the  function  Rizy  AO  is  the  autocorrelation 
function  of  the  instantaneous  mutual  coherence  of  the  two 
waves,  i.e., 

R(z,  A t)  =  (F(z,  t  +  A t)F*(z,  0>.  (2D 

As  the  time  delay  A  t  approaches  infinity,  the  two  random 
variables  F(z ,  t  +  A*)  and  F(z,  0  become  uncorrelated. 
According  to  Eq.  (21),  we  have 

R(z,  oo)  =  \(F(Zy  t)) i2.  (22) 

In  other  words,  R(z ,  oo)  is  the  squared  mutual  coherence 
of  the  two  waves.  Using  Eq.  (18),  we  can  express  R(zy  °0 
explicitly. 

With  Eq.  (20),  the  problem  of  solving  for  the  signal  in¬ 
tensity  is  reduced  to  the  problem  of  solving  for  the  auto¬ 
correlation  function  of  the  instantaneous  mutual  coher¬ 
ence  of  the  two  waves.  Using  Eqs.  (17)  and  (21),  we 
obtain 


dR(z,  M) 


=  ri  f+0 
2  r  J.. 


R(z ,  A t  +  AO)expl - —  IdAO. 


R(0y  AO  -  /1(0)/2(0)(exp[i Bit  +  A*  +  td)  -  i&it  +  AO 
-  i6(t  +  td)  +  iO(t)]).  (26) 

The  above  expression  involves  the  random  phase  0(0  at 
four  different  times:  t  +  Af  +  td ,  t  +  A tyt  +  ^  ,  and  f . 
The  temporal  sequence  of  these  times  plays  an  important 
role  in  determining  the  explicit  expression  of  the  function 
i?(0,  AO.  As  is  shown  in  Fig.  2,  there  are  four  possible 
sequences  for  td>  0  and  four  possible  sequences  for  td 
<0.  As  an  example,  we  consider  one  such  sequence,  the 
case  A t  >  td>  0.  The  random  phase  shift  A d'(td) 
=  Bit  +  A*  +  td)  -  Bit  +  AO  and  the  random  phase 
shift  A &’{td)  =  Bit  +  td)  -  Bit)  are  over  nonoverlapping 
time  intervals,  and  therefore  they  are  statistically  inde¬ 
pendent.  Thus  we  can  simplify  Eq.  (26)  as,  for  this  par¬ 
ticular  case, 

#(0,  AO  =  /i(O)/2(O){exp[iA0/(^)]){exp[-iA0,,(^)]). 

(27) 


We  define  the  Fourier  transform  of  the  function  Rizy  AO 


RiZy  &  = 


Rizy  AOexp(ffAOdAf.  (24) 


Using  Eqs.  (23)  and  (24),  we  obtain 


R(z,  £)  =  i*(0,  Oexp 


i  +  e^y 


Given  the  boundary  condition  RiOy  AO,  we  can  obtain 
Riz ,  AO  with  Eqs.  (24)  and  (25). 

Using  Eqs.  (10)  and  (11),  the  boundary  condition 
Ri 0,  AO  can  be  written  as 


Using  Eq.  (12)  and  the  assumption  that  the  random  phase 
shift  A  Bi  St)  over  an  arbitrarily  given  time  period  St  has  a 
Gaussian  distribution,  we  obtain 


R(0,t)  =  71(O)/2(O)exp[-(A0'«d)2)/2] 

X  exp[  -(A  0"(td)2)/2] 

=  Ji(0)/2(0)exp(-<Jarfd).  (28) 

The  boundary  condition  i?(0,  AO  for  the  other  cases  can 
be  obtained  similarly.  The  general  result  for  the  com¬ 
plete  boundary  condition  can  be  written  as 


At>td>0 


0>td>At 


t  t+td  t+At  t+At+td 


t+At+td  t+At  t+td  1 


td>At>0 


0>At>td 


t  t+At  t+td  t+At+td 


t+At+td  t+td  t+At  1 


td>-At>0 


0>-At>td 


t+At  1  t+At+td  t+td 


t+At+td 


-At>td>0 


0>td>-At 


t+At  t+At+td  1 


t  t+At+td  t+At 


Fig.  2.  Eight  possible  time  sequences  of  t  +  A£  +  td ,  t  +  A t,  t  +  td ,  and  t. 
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72( 0,  AO 


f/1(0)/2(0)exp(-^|^|) 

i/1(0)J2(0)exp(-^|A^) 


>  m 


By  use  of  Eqs.  (17),  (20),  (24),  (25),  and  (29),  the  signal  in¬ 
tensity  and  the  mutual  coherence  of  the  two  waves  can  in 
general  be  obtained  numerically. 

B.  Discussion 

In  this  subsection  we  discuss  some  general  characteristics 
and  some  special  cases  that  can  be  solved  analytically. 
We  first  describe  an  interesting  relationship  between  the 
signal  intensity  and  the  mutual  coherence  of  the  two 
waves. 

In  general,  the  signal  wave  and  the  pump  wave  are 
partially  coherent  at  the  incidence  plane.  We  can  write 
the  signal  wave  amplitude  at  the  incidence  plane  as  the 
sum  of  two  terms,  i.e., 

®i(0,f)  =  Elc(0,t)  +Eu(0,t),  (30) 

where  the  term  7£lc(0,  t)  is  fully  coherent  with  the  pump 
wave  at  the  incidence  plane  and  the  term  Eu( 0,  t)  is  to¬ 
tally  incoherent  with  the  pump  wave  at  the  incidence 
plane.  This  partition  is  unique.  The  two  terms  can  be 
written  as 

<F(0,O> 

Elc(0,t)  =  X^-W-S2( 0,0,  (3D 

(F( 0,  0) 

Eu( 0,  0  =  E i(0,  0  -  ^ -  E2(01  O.  (32) 

With  Eqs.  (30)-(32),  we  have  separated  the  signal  wave, 
before  it  enters  the  photorefractive  medium,  into  two  in¬ 
dependent  wave  components.  One  is  fully  coherent  with 
the  pump  wave  and  the  other  is  totally  incoherent  with 
the  pump  wave.  According  to  Eqs.  (7)-(9),  the  coupled- 
wave  equations  are  linear  equations  of  the  signal-wave 
amplitude.  Therefore  these  two  signal-wave  components 
will  also  propagate  through  the  photorefractive  medium 
coupling  with  the  pump  wave  independently.  The  ampli¬ 
tudes  of  two  signal-wave  components  throughout  the  pho¬ 
torefractive  medium  are  denoted  Elc(z ,  t)  and  Eu(z,  O, 
respectively.  The  intensities  of  the  two  signal-wave  com¬ 
ponents  are  denoted  7lc(z)  and  Iu(z),  respectively.  The 
autocorrelation  function  of  the  instantaneous  mutual  co¬ 
herence  of  the  two  signal-wave  components  with  the 
pump  wave  are  denoted  Rc(zf  A t)  and  72, (z,  AO,  respec¬ 
tively.  Using  Eqs.  (31)  and  (32),  we  obtain 

/i(0)  =  /lc(0)  +  Iu(  0),  (33) 

R( 0,  AO  =  Rc( 0,  AO  +  72*(0,  AO.  (34) 

Note  that  the  ensemble  averages  of  the  cross-product 
terms  in  Eqs.  (33)  and  (34)  are  all  zero.  According  to 
their  definitions,  72c( 0,  AO  is  a  constant  of  the  time  delay 
A 0  whereas  72j(0,  AO  approaches  zero  as  the  time  delay 
A t  goes  to  infinity.  Thus,  given  the  function  72(0,  AO,  we 
can  obtain  72c(0,  AO  and  72;(0,  AO  as 

72c(0,  AO  =  72(0,oo),  (35) 

72,(0,  AO  =  72(0,  AO  -  72(0,  oo).  (36) 


The  propagation  of  the  signal  wave  component  that  is 
fully  coherent  with  the  pump  wave  at  the  incidence  plane 
can  be  solved  explicitly.  The  boundary  condition  is  given 
by  72c(0,  AO  =  7Zc(0)72(0).  It  is  a  constant  of  the  time 
delay  AO  and  its  Fourier  transform  is  a  delta  function. 
Using  Eqs.  (24)  and  (25),  we  obtain 

72c(z,  AO  =  7lc(0)/2(0)exp(yz).  (37) 

Note  that  the  function  72c(z,  AO  is  a  constant  of  the  time 
delay  A*  at  any  position  z.  The  signal- wave  intensity 
can  thus  be  expressed  as 

7lc(z)  =  7lc(0)exp(yz).  (38) 

With  Eqs.  (37)  and  (38),  we  can  prove  that  the  absolute 
value  of  the  normalized  mutual  coherence  of  the  two 
waves  is  1  throughout  the  photorefractive  medium.  The 
signal- wave  amplitude  can  be  written  as  Elc(z,  O 
=  7£lc(0,  t  -  z/u)exp(yz/2).  The  signal  intensity  gain  is 
the  same  as  in  the  case  in  which  the  two  waves  are  mono¬ 
chromatic  waves.  The  normalized  spectrum  of  the  signal 
wave  does  not  change  in  the  two-wave  mixing  process  and 
remains  the  same  as  the  pump  wave.  Note  that  the 
propagation  of  the  signal-wave  component  that  is  coher¬ 
ent  with  the  pump  wave  is  independent  of  the  photore¬ 
fractive  relaxation  time  r  when  the  steady  state  is 
reached. 

The  above  analysis  shows  that  the  signal-wave  compo¬ 
nent  that  is  totally  coherent  with  the  pump  wave  at  the 
incidence  plane  remains  totally  coherent  with  the  pump 
wave  throughout  the  photorefractive  medium.  In  addi¬ 
tion,  according  to  Eq.  (18),  the  signal-wave  component 
that  is  totally  incoherent  with  the  pump  wave  at  the  in¬ 
cidence  plane  remains  totally  incoherent  with  the  pump 
wave  throughout  the  photorefractive  medium.  There¬ 
fore,  Eqs.  (31M36)  apply  not  only  at  the  incidence  plane 
but  also  throughout  the  photorefractive  medium. 

Since  the  propagations  of  the  two  signal-wave  compo¬ 
nents  are  independent  of  each  other  and  since  the  propa¬ 
gation  of  the  signal-wave  component  that  is  totally  coher¬ 
ent  with  the  pump  wave  can  be  easily  determined,  we 
now  consider  the  propagation  of  the  signal  wave  compo¬ 
nent  that  is  totally  incoherent  with  the  pump  wave. 

Without  loss  of  generality,  we  consider  the  boundary 
condition  given  by  Eq.  (29).  We  assume  that  the  time  de¬ 
lay  td  between  the  signal  wave  and  the  pump  wave  is 
much  larger  than  the  coherence  length  of  the  source  laser 
wave  so  that  the  two  waves  are  totally  incoherent.  Un¬ 
der  this  assumption,  the  boundary  condition  [Eq.  (29)] 
can  be  simplified  as 

72(0,  AO  =  71(0)72(0)exp(-<MA*|),  (39) 

and  its  Fourier  transform  is 


72(0,  O  =  7,(0172(0)  +-  j2.  (40) 

The  function  72(0,  f)  has  a  Lorentzian  shape.  Its  maxi¬ 
mum  is  72(0,  0)  =  Il{Q)I2(0)/(7rS(jj).  Its  FWHM  is  2<Sw, 
i.e.,  twice  the  bandwidth  of  the  source  laser  wave.  Gen¬ 
erally  speaking,  when  the  signal  wave  and  the  pump 
wave  are  incoherent  the  boundary  condition  72(0,  f)  is  a 
function  of  £  with  finite  height  and  finite  width;  when  the 


(36) 
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signal  wave  and  the  pump  wave  are  coherent,  the  bound¬ 
ary  condition  R( 0,  £)  is  a  delta  function  of  $ 

We  now  describe  some  general  characteristics  of  the 
case  in  which  the  two  waves  are  totally  incoherent.  Look 
at  the  function  exp[yz/(l  +  £2t2)]  in  Eq.  (25).  It  is  an  in¬ 
creasing  function  of  yz.  According  to  Eqs.  (24)  and  (25), 
both  the  function  R(z ,  £)  and  the  function  R(z ,  A t)  are 
increasing  functions  of  yz.  Thus,  according  to  Eq.  (20), 
the  signal  intensity  will  be  amplified  even  if  the  two 
waves  are  incoherent.  The  function  exp[yz/(l  +  f  2r2)]  is 
a  decreasing  function  of  the  photorefractive  relaxation 
time  r.  According  to  Eqs.  (20),  (24),  and  (25),  the  signal 
intensity  gain  decreases  as  r  increases.  The  function 
exp[yz/(l  +  £2t2)]  is  also  a  decreasing  function  of  |£|.  Its 
FWHM  with  respect  to  the  argument  £  is  a  decreasing 
function  of  yz.  When  yz  >  In  2,  the  width  is  approxi¬ 
mately  [In  2 l{yz)]mlr.  According  to  Eqs.  (24)  and  (25), 
the  FWHM  of  the  function  R(z ,  A*)  with  respect  to  the 
argument  A t  increases  as  yz  increases.  In  other  words, 
the  correlation  time  of  the  instantaneous  mutual  coher¬ 
ence  of  the  two  waves  increases  as  the  signal  wave  be¬ 
comes  amplified. 

We  further  consider  some  special  cases  that  can  be 
solved  analytically. 

First,  we  consider  the  case  in  which  the  photorefractive 
grating  relaxation  time  constant  is  much  shorter  than  the 
coherence  time  of  the  two  interacting  waves,  _i.e.,  Sur 
<  1.  In  this  case,  the  width  of  the  function  R( 0,  f)  is 
much  smaller  than  the  width  of  the  function  exp[yz/(l 
4-  £2t2)],  i.e.,  6<dt<  [in  2/(4 yz)]172.  Note  that 

[In  2/(4  yz)] 1/2  is  of  the  order  of  unity  even  for  a  relatively 
large  value  of  yz.  Therefore,  the  function  exp[yz/(l 
+  £2t2)]  in  Eqs.  (24)  and  (25)  can  be  replaced  with 
exp(yz),  i.e.,  r  =  0.  This  leads  to 

R(z ,  M)  =  R( 0,  Af)exp(yz),  (41) 

.  h(z)  =  /i(0)exp(yz).  (42) 

The  signal  intensity  gain  is  the  same  as  in  the  case  of  the 
two  waves’  being  totally  coherent.  The  function  R (z ,  At) 
retains  its  temporal  profile  while  being  amplified.  In  this 
case  the  photorefractive  grating  follows  the  fluctuation  of 
the  partially  coherent  optical  waves  instantaneously,  i.e., 
Q(z,  t)  =  j Ex(z,  t')E2*(z ,  t’)  according  to  Eq.  (7),  and 
the  signal  wave  amplitude  can  be  directly  obtained  from 
Eq.  (8)  as  Ex(z,  t)  =  JSjfO,  t  -  z/v)exp(yz/2).  Therefore 
the  normalized  spectrum  of  the  signal  wave  does  not 
change  in  the  two-wave  mixing  process  and  does  not  have 
to  be  the  same  as  the  pump  wave.  This  is  in  agreement 
with  the  fast  response  of  the  material.  Since  the  signal 
intensity  gain  decreases  as  the  photorefractive  relaxation 
time  r  increases  when  the  two  waves  are  totally  incoher¬ 
ent,  the  signal  intensity  gain  when  the  two  waves  are  to¬ 
tally  incoherent  will  in  general  be  smaller  than  when  the 
two  waves  are  totally  coherent  and  the  photorefractive  re¬ 
laxation  time  r  is  not  negligible. 

Second,  we  consider  the  case  in  which  the  photorefrac¬ 
tive  grating  relaxation  time  constant  is  much  longer  than 
the  coherence  time  of  the  two  interacting  waves,  i.e., 
Swt  >  1.  In  this  case  the  width  of  the  function  R{  0,  £)  is 
much  larger  than  the  width  of  the  function  exp[yz/(l 
4-  £2t2)].  As  a  result,  the  boundary  condition  R(0>  £)  in 
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Eq.  (25)  can  be  replaced  with  f?(0,  0)  =  /1(0)/2(0)/ 
(ttSoj).  When  yz  <  1,  the  function  exp[yz/(l  +  £2t2)] 
can  be  replaced  with  1  +  yz/(  1  +  £2t2).  This  leads  to 

yz  I  \At\\ 

R(z,  M)  =  R(0,  At)  +  —  exp  | — —  j,  (43) 

/,(.)- /,(0)(l+  £;).  <44> 

The  signal  intensity  gain  coefficient  is  yzl{Scor)  in  the 
low-gain  limit.  The  incremental  change  of  the  function 
R(zy  A t),  i.e.,  the  second  term  on  the  right-hand  side  of 
Eq.  (43),  has  a  FWHM  of  r  In  2,  which  is  much  larger 
than  the  FWHM  of  the  function  tf(0,  At)y  i.e.,  In  2 ISo). 
When  yz  >  1,  the  function  exp[yz/(l  +  f2r2)]  can  be  re¬ 
placed  with  exp[yz(l  -  f2r2)].  This  leads  to 


1  exp  (yz)  j 

M2\ 

(45) 

R(z,  M)  = 

S»r  ^  GXP1 

4  yzT2/ 

^i(0)  exp('yz) 

(46) 

Hz)  = 

dtoT  yj  Tjyz 

In  the  high-gain  limit,  the  signal  intensity  gain  coefficient 
is  close  to  yz,  the  same  as  in  the  case  in  which  the  two 
waves  are  totally  coherent,  whereas  the  signal  intensity 
gain  is  at  least  Sojt  times  smaller  than  in  the  case  in 
which  the  two  waves  are  totally  coherent.  The  temporal 
profile  of  the  function  R(z,  At)  is  Gaussian  with  a  FWHM 
of  2  rypyz.  Codirectional  two-wave  mixing  with  partially 
coherent  waves  was  studied  by  Bogodaev  et  al.  under  the 
assumption  that  the  amplitude  of  the  photorefractive 
grating  is  temporally  stationary  and  proportional  to  the 
mutual  coherence  of  the  two  interacting  waves,  i.e., 
Q(z,  t)  =  (F(z,  £)),  in  the  case  Scot>  l.7  Under  the 
undepleted  pump  approximation,  their  results  show  that 
the  signal  wave  does  not  become  amplified  at  all  when  it 
is  incoherent  with  the  pump  wave,  whereas  the  signal  in¬ 
tensity  and  the  mutual  coherence  of  the  two  waves  are 
the  same  as  those  given  by  Eqs.  (18)  and  (38)  when  the 
two  waves  are  totally  coherent.  The  assumption  of  tem¬ 
porally  stationary  photorefractive  grating  is  valid  when 
the  amplification  of  the  signal  intensity  due  to  a  signal- 
wave  component  incoherent  with  the  pump  wave  is  negli¬ 
gible  compared  with  that  due  to  a  signal-wave  component 
that  is  coherent  with  the  pump  wave. 

C.  Simulations 

We  now  show  some  numerical  results  for  the  signal  inten¬ 
sity  gain  using  the  boundary  condition  given  by  Eq.  (29). 
In  Fig.  3,  we  show  the  signal  intensity  gain  as  a  function 
of  the  position  z  in  the  photorefractive  medium  for  a  given 
time  delay  td  between  the  signal  wave  and  the  pump 
wave  at  the  incidence  plane  and  for  five  different  values  of 
the  photorefractive  grating  relaxation  time.  The  simula¬ 
tion  parameters  are  interaction  length  L  =  1.0  cm,  index 
of  refraction  n  =  2.3,  photorefractive  coupling  coefficient 
y  =  1.0  cm-1,  source  laser  wave  bandwidth  Sw 
=  27rXl.O  X  106  s“\  time  delay  td  -  0.3  /xs,  and  five 
different  values  of  the  photorefractive  grating  relaxation 
time  r  =  0.0,  0.4,  1.0,  2.5,  and  «>  M s .  For  the  three  finite 
cases,  Swt  is  of  the  order  of  unity.  According  to  Fig.  3, 


X.  Yi  and  P.  Yeh 

the  signal  intensity  gain  increases  as  functions  of  2  in  the 
photorefractive  medium  for  a  given  photorefractive  relax¬ 
ation  time  r  and  increases  as  the  photorefractive  relax¬ 
ation  time  r  decreases.  The  signal  intensity  is  the  sum  of 
the  intensities  of  the  signal-wave  components  that  are  co¬ 
herent  with  the  pump  wave  and  of  those  incoherent  with 
the  pump  wave.  The  intensity  gain  of  the  signal-wave 
component  that  is  coherent  with  the  pump  wave  is  inde¬ 
pendent  of  the  photorefractive  relaxation  time  r.  There¬ 
fore,  according  to  Fig.  3,  the  intensity  gain  of  the  signal- 
wave  component  that  is  incoherent  with  the  pump  wave 
increases  as  the  photorefractive  relaxation  time  r  de¬ 
creases.  When  the  photorefractive  relaxation  time  r  is 
zero,  the  intensity  gains  of  the  two  signal-wave  compo¬ 
nents  are  equal.  When  the  photorefractive  relaxation 
time  r  is  not  zero,  the  intensity  gain  of  the  signal  wave 
component  that  is  incoherent  with  the  pump  wave  is 
smaller  than  that  of  the  signal  wave  component  that  is  co¬ 
herent  with  the  pump  wave. 

In  Fig.  4  we  show  the  signal  intensity  gain  at  the  exit 
plane  as  a  function  of  the  time  delay  td  between  the  two 
waves  at  the  incidence  plane  for  five  different  values  of 
the  photorefractive  grating  relaxation  time.  The  simula¬ 
tion  parameters  are  the  same  as  those  of  Fig.  3.  Since 
the  curves  of  the  signal  intensity  gain  versus  the  time  de¬ 
lay  td  are  symmetric  with  respect  to  the  axis  when  the 
time  delay  td  is  zero,  we  show  only  the  curves  for  the  posi¬ 
tive  time  delay  td .  We  notice  that  the  signal  intensity 
gains  are  decreasing  functions  of  the  time  delay  td .  For 
all  cases,  the  signal  intensity  gains  reach  asymptotic  val¬ 
ues  when  the  time  delay  td  is  a  couple  of  the  coherence 
time  of  the  source  laser  wave  or  larger.  When  Scot  is 
much  larger  than  1,  there  is  no  coupling  between  the  two 
waves  if  the  two  waves  are  incoherent  at  the  incidence 
plane  owing  to  a  large  time  delay  (e.g.,  td  >  1  /as  in  Fig. 
4).  When  Scot  is  of  the  order  of  unity  or  less,  there  is 
always  coupling  between  the  two  waves,  even  if  the  two 
waves  are  completely  incoherent  at  the  incidence  plane 
owing  to  a  large  time  delay  (e.g.,  td  >  1  fiis  in  Fig.  4). 
Also  according  to  Fig.  4,  the  signal  intensity  gain  in¬ 
creases  as  the  photorefractive  grating  relaxation  time  de¬ 
creases.  The  increase  of  the  signal  intensity  gain  with 
the  decrease  of  the  photorefractive  grating  relaxation 


Fig.  3.  Signal  intensity  gain  as  a  function  of  the  position  in  the 
photorefractive  medium  for  a  given  time  delay  td  between  the 
signal  wave  and  the  pump  wave  at  the  incidence  plane  and  for 
five  different  values  of  the  photorefractive  grating  relaxation 
time. 
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Fig.  4.  Signal  intensity  gain  at  the  exit  plane  as  a  function  of 
the  time  delay  td  between  the  two  waves  at  the  incidence  plane 
for  five  different  values  of  the  photorefractive  grating  relaxation 
time. 

time  becomes  more  significant  when  the  time  delay  td  is  a 
couple  of  the  coherence  time  of  the  source  laser  wave  or 
larger.  The  dependence  of  the  signal  intensity  gain  on 
the  time  delay  td  between  the  two  waves  at  the  incidence 
plane  can  be  easily  understood  as  follows.  Increasing  the 
time  delay  td  between  the  two  waves  increases  the  rela¬ 
tive  weight  of  the  signal-wave  component,  which  is  inco¬ 
herent  with  the  pump  wave.  This  explains  the  decrease 
of  the  signal  intensity  gain.  When  the  time  delay  td  is  a 
couple  of  the  coherence  time  of  the  source  laser  wave,  the 
signal  wave  becomes  completely  incoherent  with  the 
pump  wave  and  the  signal  intensity  gain  reaches  its 
asymptotic  value. 

In  this  section  the  stochastic  coupled-wave  equations 
are  reduced  to  nonstochatic  differential  equations,  from 
which  the  signal  intensity  and  the  mutual  coherence  of 
the  two  waves  can  be  obtained.  The  autocorrelation 
function  of  the  instantaneous  mutual  coherence  of  the  two 
waves  at  the  incidence  plane  is  found  to  be  the  necessary 
and  sufficient  boundary  condition  to  determine  the  evolu¬ 
tion  of  the  signal  intensity  through  the  two-wave  mixing 
process.  Our  results  indicate  that  both  the  signal-wave 
component  that  is  coherent  with  the  pump  wave  and  the 
signal-wave  component  that  is  incoherent  with  the  pump 
can  become  amplified  while  maintaining  their  respective 
mutual  coherence  properties  with  the  pump  wave.  The 
intensity  gain  of  the  signal-wave  component  that  is  inco¬ 
herent  with  the  pump  wave  depends  on  the  photorefrac¬ 
tive  relaxation  time  constant,  whereas  the  signal-wave 
component  that  is  coherent  with  the  pump  wave  does  not. 
The  effect  of  linear  absorption  can  be  easily  taken  into  ac¬ 
count. 


4.  SIGNAL  SPECTRUM 

In  this  section  we  solve  for  the  signal  spectrum.  The  sto¬ 
chastic  coupled  wave  equations  are  rewritten  into  a  series 
form,  from  which  an  asymptotic  solution  of  the  signal- 
wave  spectrum  can  be  obtained. 

Since  the  signal- wave  amplitude  E i(z,  t)  is  a  station¬ 
ary  random  process,  the  signal  spectrum  denoted 
En(z ,  00)  can  be  mathematically  defined  as 
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1  f +x 

En(z,  a>)  =  (E1(z,  t)Ei*(z,  t  +  Af)) 

/7T  J  _o o 

X  exp[-i(o>  ~  o>0)A£]dA£,  (47) 

where  {E^z,  t)Ex*{z ,  t  +  AO)  as  a  function  of  A£  is  the 
autocorrelation  function  of  the  signal  wave  amplitude 
Ei(z,  t).  According  to  the  Wiener-Khintchine  theorem, 
En(z ,  co)do)  is  the  average  signal  intensity  that  is  due  to 
frequency  components  between  (d  and  w  +  dw.  Since 
the  signal-wave  spectrum  and  the  signal-wave  autocorre¬ 
lation  function  are  Fourier  transform  pairs,  we  can  obtain 
the  signal  spectrum  by  solving  for  the  signal-wave  auto¬ 
correlation  function  ( Ei(z ,  t)Ei*(z ,  t  +  AO). 

The  coupled-wave  equations  (7)-(9)  can  be  expanded 
into  a  series  of  equations.  First,  we  write  the  signal- 
wave  amplitude  Ex{z,  t)  into  the  summation  of  a  series 
as 

Ex(z ,  t)  =  E10(z ,  t)  +  En(z ,  t)  +  E12(zy  t)  +  ...  . 

(48) 

Each  summation  term  in  Eq.  (48)  is  determined  by  the 
boundary  condition  and  the  series  expansion  of  the 
coupled  wave  equations  (7)-(9)  as  will  be  shown  by  Eqs. 
(50)-(54).  We  then  write  the  grating  amplitude  Q(z,  t) 
into  the  summation  of  a  series  as 

Q(z,  t)  =  Qio(2>  0  Qu(zt  0  4-  Q  12(2 1  0  +  •••  > 

(49) 

where  Qy(z ,  0  is  the  component  due  toE^z,  0  and  the 
pump  wave  E2(zy  0.  We  do  not  expand  the  pump-wave 
amplitude  E2(zf  0,  which  is  a  constant  in  our  nonde- 
pleted  pump  approximation.  Each  component  of  the 
grating  amplitude  in  Eq.  (49)  is  determined  by 


dExi(z,  t)  1  dEv(z,  t)  y  Quj-iM,  t) 

Tz  +  Z  Jt  =  2  7^(oj  £z(  ’ f)' 

(54) 

By  summing  both  sides  of  Eq.  (54)  over  j  —  0, 1,  2, we 
obtain  Eq.  (8).  Therefore  Eqs.  (48)-(54)  are  the  series  ex¬ 
pansion  form  of  the  coupled-wave  equations. 

Using  Eq.  (48),  we  can  expand  the  signal-wave  autocor¬ 
relation  function  (E^z,  t)E{*(z,  t  +  &))  into  a  series 
summation  as 

(Ex(z,  t)E^(Zy  t  +  AO) 

-  (El0(z,  t)E10*(zf  t  +  AO) 

+  [(El0(zy  t)En*(z,  t  +  AO) 

+  {Eu(zy  t)E10*(z,  t  +  AO)]  +  •  (55) 

For  example,  the  single  product  term  [E  im 
x  (z,  t)Ein*(z,  t  +  AO]  belongs  to  the  (m  +  n)th  order 
term  of  the  autocorrelation  function  (Ex{zy  0 
x  Ei*(zy  t  +  AO).  The  product  terms  in  Eq.  (55)  can 
be  regrouped  so  that  the  autocorrelation  function 
(Ei(Zy  t)Ei*(zt  t  +  AO)  is  expanded  into  a  Taylor’s  se¬ 
ries  of  yz. 

The  zeroth-order  term  of  the  Taylor  series  is 

{E10(z,  t)El0*(z,  t  +  At))  =  /].(0)exp|  ~  Aw|At|). 

(56) 

Equation  (56)  reflects  that  the  source  laser  wave  has  been 
assumed  to  have  a  Lorentzian  line  shape.  The  first-order 
term  of  the  Taylor  series  contains  two  single-product 
terms.  These  terms  can  be  obtained  from 


r  dQlj(*’  -  +  Quiz,  t)  =  Eu(z,  t)E2*(z,  t),  (50) 

where  the  index  j  =  0, 1,  2, ... .  With  each  of  the  compo¬ 
nents  determined  this  way,  the  grating  amplitude  will  au¬ 
tomatically  satisfy  Eq.  (6)  or  Eq.  (7).  Physically,  this 
means  that  the  j'th-order  term  of  the  grating  amplitude  is 
generated  by  the  jth-order  term  of  the  signal  wave  and 
the  pump  wave.  The  boundary  conditions  at  the  inci¬ 
dence  plane  for  the  components  of  the  signal-wave  ampli¬ 
tude  in  Eq.  (48)  are  given  by 

£io(0,  t)  ~  Ei(0y  t)y  (51) 

E.jiOyt)  =  0,  (52) 

where  j  =  1,  2, ...  .  In  this  way  the  boundary  conditions 
at  the  incidence  plane  for  the  signal-wave  amplitude  are 
also  automatically  satisfied.  The  zeroth-order  term  of 
the  signal  wave  is  assumed  to  pass  through  the  photore- 
fractive  medium  without  being  affected  by  the  photore- 
fractive  grating,  i.e., 

dE iq(z y  t)  1  SE10(Zy  t) 

-  +  77 - ^ -  =  0-  (5o) 

dz  V  dt 

For  j  >  0,  the  j th  order  term  of  the  signal  wave  is  as¬ 
sumed  to  be  generated  by  the  pump  wave  reflection  from 
the  (j  -  l)th  order  term  of  the  photorefr active  grating, 
i.e., 


(Eu(x,  t)E io* (z,  t  +  A t)) 

-  T  m~r  />■<<.•  <  +  < +  4,,) 

X  E2(z,  t)Eio*(z,  t  +  Af))exp|—  JdAf',  (57) 

(Ew(z,  t)En*(z,  t  +  AO) 

=  (Eu{z,  t)E10*(z,  t  -  AO)*.  (58) 

Using  the  boundary  condition  given  by  Eq.  (52),  we  can 
show  that  both  sides  of  Eq.  (57)  are  zero  at  the  incidence 
plane  z  =  0.  Using  the  series  expansion  form  of  the 
coupled-wave  equations,  we  can  show  that  the  first  de¬ 
rivatives  over  the  position  z  of  both  sides  of  Eq.  (57)  are 
equal.  This  proves  the  validity  of  Eq.  (57).  According  to 
Eq.  (57),  the  coefficient  of  the  first-order  term  of  the  Tay¬ 
lor  series  is  related  to  a  one-layer  integration  of  a  fourth- 
order  coherence  function  of  the  two  wave  amplitudes  at 
the  incidence  plane.  When  the  source  laser  wave  is  as¬ 
sumed  to  have  only  phase  noise,  the  fourth-order  coher¬ 
ence  function  in  Eq.  (57)  can  be  obtained  in  the  same  way 
as  we  obtain  the  fourth-order  coherence  function  i2(0,  A t) 
given  in  Eq.  (29).  Its  general  expression  is  quite  compli¬ 
cated.  When  the  time  delay  td  between  the  signal  wave 
and  the  pump  wave  is  further  assumed  to  be  much  larger 
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than  the  coherence  time  of  the  source  laser  wave,  the 
fourth-order  coherence  function  in  Eq.  (57)  can  be  ex¬ 
pressed  as 

(E10(z,  t  +  Af ')E2*(z,  t  +  A t')E2(z,  t)E10*(z,  t  +  Af)) 


=  /1(0)/2(0)exp 


--  -  Af|  +  | A£ ' |) 


(59) 


Using  Eq.  (59),  we  obtain 


(E10{z,  t)En*(z,  t  +  At))  +  (E n(z,  O£}0*(z,  t  +  Af)) 


=  /j(0)exp  --  5&j|At| 


yz 


SCJT  +  1 


X 


1+- 


S(i)T 


(60) 


Fig.  5.  First  two  orders  of  the  Taylor’s  series  of  the  autocorre¬ 
lation  function  (Ex{z,  t)E i*U,  t  +  A*))  as  a  function  of  the  time 
delay  M. 


In  the  case  Scjt  <  1,  Eq.  (60)  can  be  simplified  as 
<Eio(z,  t)En*(z,  t  +  At))  +  <£n(z,  t)Ew*(z,  t  +  At)) 

=  71(0)exp|  —  —  <Sco|Af|j  yz.  (61) 

In  the  case  S(dt>  1,  Eq.  (60)  can  be  simplified  as 
(E10(z,  t)En*(z,  t  +  At))  +  (En(z,  t)E  10* (z,  t  +  At)) 

=  7i(0)exp|  —  5<t)|Af||  (62) 

According  to  Eqs.  (56),  (61),  and  (62),  when  the  signal 
wave  and  the  pump  wave  are  independent  at  the  inci¬ 
dence  plane  and  the  gain-length  product  yz  is  small  com¬ 
pared  with  unity,  the  normalized  spectrum  of  the  signal 
wave  does  not  change  in  the  cases  1  and  Swt 

>  1. 

In  general,  there  are  j  +  1  single  product  terms  in  the 
j th-order  term  of  the  Taylor  series.  The  coefficient  of  the 
j th-order  term  of  the  Taylor  series  is  related  to  j  -level  in¬ 
tegrations  of  (j  +  l)th-order  coherence  functions  of  the 
two  wave  amplitudes  at  the  incidence  plane.  The  coher¬ 
ence  functions  of  the  two  wave  amplitudes  at  the  inci¬ 
dence  plane  can  be  obtained  from  the  boundary  conditions 
in  the  same  way  as  we  obtain  the  fourth-order  coherence 
function  #(0,  M)  given  in  Eq.  (29).  As  we  have  shown, 
the  signal-wave  intensity  increases  as  exp(>e).  Since  the 
signal-wave  intensity  is  equal  to  the  autocorrelation  func¬ 
tion  (Ex(z,  t)Ex*(z ,  t  +  AO)  at  the  time  delay  A t  =  0, 
we  expect  the  correlation  function  to  increase  as  exp(yz) 
for  the  time  delay  A£  ^  0  as  well.  Therefore  the  number 
of  orders  in  the  Taylor’s  series  that  is  needed  to  achieve  a 
certain  accuracy  increases  quickly  as  yz  increases.  The 
series  expansion  method  works  most  efficiently  when  yz 
is  close  to  1  or  smaller.  Fortunately,  yz  is  usually  small 
for  high-speed  photorefractive  materials. 

For  example,  when  yz  is  small  compared  with  unity, 
only  the  first  two  orders  of  the  Taylor’s  series  need  to  be 
taken  into  account.  In  Fig.  5  we  show  the  first  two  orders 
of  the  Taylor’s  series  separately  as  a  function  of  the  time 
delay  A t.  Both  functions  are  normalized  by  their  respec¬ 
tive  maximum  values.  The  time  delay  between  the 
two  waves  at  the  incidence  plane  is  assumed  to  be  much 
larger  than  the  coherence  time  of  the  source  laser  wave. 
In  this  particular  case,  according  to  Eqs.  (55)  and  (59),  the 


maximum  value  of  the  zero  th-order  term  is  I\(0)f  and  the 
maximum  value  of  the  first-order  term  is  yzIx(0)/(So)T 
+  1).  The  bandwidth  of  the  source  laser  wave  is  as¬ 
sumed  to  be  1  MHz.  The  photorefractive  relaxation  time 
is  assumed  to  be  1  /is.  According  to  Fig.  5,  the  widths  of 
the  first  two  orders  of  the  autocorrelation  function  are  not 
equal.  The  bandwidth  of  the  signal  wave  will  become 
narrower  as  the  signal  wave  propagates  through  the  pho¬ 
torefractive  medium. 

In  this  section  we  have  presented  an  asymptotic  solu¬ 
tion  of  the  signal  spectrum.  The  solution  is  relatively 
simple  in  the  low-gain  limit.  When  the  coherence  time  of 
the  source  laser  wave  is  comparable  with  the  photorefrac¬ 
tive  relaxation  time,  the  signal  wave  spectrum  will  be¬ 
come  narrower  as  the  signal  wave  becomes  amplified. 
When  the  coherence  time  of  the  source  laser  wave  is  much 
longer  or  much  shorter  than  the  photorefractive  relax¬ 
ation  time,  the  signal  wave  spectrum  remains  unchanged. 

5.  CONCLUSIONS 

We  have  investigated  codirectional  two-wave  mixing  with 
partially  coherent  waves  in  high-speed  photorefractive 
media.  When  the  spectral  bandwidth  of  the  optical 
waves  is  close  to  the  speed  of  the  photorefractive  media, 
i.e.,  Scot  is  of  the  order  of  unity,  the  photoinduced  index 
grating  will  fluctuate  with  the  amplitudes  of  the  optical 
waves.  The  stochastic  coupled-wave  equations  are 
solved  in  the  nondepleted  pump  region  with  the  assump¬ 
tion  that  the  pump  wave  has  only  phase  noise.  Results 
for  the  evolution  of  the  signal  intensity,  the  mutual  coher¬ 
ence  of  the  two  waves,  and  the  signal  spectrum  are  ob¬ 
tained.  We  found  that  the  two-wave  mixing  process  in 
high-speed  photorefractive  media  is  significantly  different 
from  that  in  low-speed  photorefractive  media.  In  the 
low-speed  case  we  need  to  know  only  the  intensities  and 
the  mutual  coherence  of  the  two  waves  at  the  incidence 
plane  as  the  boundary  condition  to  determine  the  evolu¬ 
tion  of  the  signal  wave  intensity,  whereas  in  the  high¬ 
speed  case  we  need  to  know  the  autocorrelation  function 
of  the  instantaneous  mutual  coherence  of  the  two  waves 
at  the  incidence  plane.  In  the  low-speed  case  a  partial 
coherence  between  the  two  waves  at  the  incidence  plane 
is  necessary  for  two-wave  mixing  and  signal-wave  ampli¬ 
fication  to  occur,  whereas  in  the  high-speed  case  two-wave 
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mixing  and  signal-wave  amplification  will  occur  even  if 
the  two  optical  waves  are  incoherent  at  the  incidence 
plane.  In  the  low-speed  case  the  two-wave  mixing  is  in¬ 
dependent  of  the  photorefractive  relaxation  time  con¬ 
stant,  whereas,  in  the  high-speed  case  the  signal  intensity 
gain  increases  as  the  photorefractive  grating  relaxation 
time  decreases  when  the  two  optical  waves  are  not  fully 
coherent  at  the  incidence  plane.  The  evolution  of  the  sig¬ 
nal  spectrum  is  also  much  more  complicated  in  the  high¬ 
speed  case  than  in  the  low-speed  case.  In  general,  the 
signal-wave  bandwidth  changes  as  a  result  of  the  two- 
waving  process.  The  analysis  in  the  high-speed  case  also 
provides  a  better  understanding  of  the  limitation  of  the 
assumptions  previously  made  in  analyzing  the  low-speed 
case. 

In  our  analysis  the  photorefractive  medium  has  been 
assumed  to  be  of  purely  diffusive  type,  which  means  that 
the  coupling  constant  y  is  a  real  number.  With  the  pres¬ 
ence  of  the  photovoltaic  effect  or  external  electric  field  in 
the  photorefractive  medium,  the  coupling  constant  y  will 
be  a  complex  number.  The  latter  case  can  be  analyzed 
similarly  to  the  former.  But  the  mathematical  treat¬ 
ments  involved  will  be  slightly  more  complicated.  The 
method  used  in  our  analysis  can  also  be  used  to  study 
stimulated  Raman  scattering  with  broadband  light. 
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We  investigate  contradirectional  two-wave  mixing  with  partially  coherent  waves  in  photorefractive  crystals. 
By  use  of  a  statistical  theory  on  linear  systems,  a  general  formulation  of  the  problem  in  the  space  and  fre¬ 
quency  domain  is  derived  and  implemented  numerically.  We  obtain  results  on  beam  intensity  and  mutual 
coherence.  The  results  on  the  enhancement  of  mutual  coherence  are  compared  with  previous  theoretical  re¬ 
sults  on  simpler  cases  and  with  experimental  measurements.  Excellent  agreements  are  achieved.  The  re¬ 
sults  also  indicate  that  the  effective  interaction  length  can  be  significantly  longer  than  the  coherence  length  of 
the  waves.  ©  1997  Optical  Society  of  America  [S0740-3224(97)01906-l] 


1.  INTRODUCTION 

Two-wave  mixing  in  photorefractive  crystals  has  been  in¬ 
vestigated  extensively  for  many  applications  including 
image  amplification,  laser-beam  cleanup,  spatial  light 
modulators,  thresholding,  and  power-limiting  devices.1,2 
Most  of  the  theoretical  works  in  this  area  are  based  on 
wave-mixing  with  mutually  coherent  waves.1,2  However, 
in  some  applications,  such  as  self-pumped  and  mutually 
pumped  phase-conjugate  mirrors3"6  and  photorefractive 
filters,7  the  effect  of  partial  temporal  coherence  in  a  two- 
wave-mixing  process  cannot  be  ignored.  Two-wave  mix¬ 
ing  with  partially  coherent  waves  has  been  studied  previ¬ 
ously  for  the  case  of  transmission-grating  interaction.8 
In  the  case  of  transmission-grating  interaction  the  optical 
path  difference  between  the  two  interacting  waves  re¬ 
mains  approximately  the  same  as  the  two  waves  propa¬ 
gating  codirectionally  through  the  photorefractive  me¬ 
dium,  especially  when  the  incident  angles  of  the  two 
waves  are  close  to  each  other.9  In  the  case  of  reflection¬ 
grating  interaction  the  optical  path  difference  between 
the  two  interacting  waves  varies  significantly  as  the  two 
waves  propagate  contradirectionally  through  the  photore¬ 
fractive  medium.  Thus  for  the  case  of  transmission¬ 
grating  interaction,  only  one  free  variable  for  the  position 
is  needed  to  describe  in  a  self-consistent  way  the  second- 
order  statistical  properties  of  the  two  optical  waves,  i.e., 
their  intensities  and  mutual  coherence,  while  at  least  two 
free  variables,  one  for  the  position  and  one  for  the  optical 
path  difference,  will  be  needed  for  the  case  of  reflection¬ 
grating  interaction.  Another  difficulty  in  studying  the 
reflection-grating  interaction  of  partially  coherent  waves 
is  to  find  a  way  to  incorporate  the  complete  boundary  con¬ 
ditions  into  the  theoretical  formulation  as  a  result  of  the 
two-point  boundary-value  problem.  In  such  a  problem  a 
complete  set  of  boundary  conditions  includes  both  the 
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second-order  self-statistical  properties  (e.g.,  self¬ 
coherence)  of  each  wave  at  its  entrance  boundary  and  the 
second-order  mutual  statistical  properties  (e.g.,  mutual 
coherence)  of  the  two  waves  at  their  respective  entrance 
boundaries.  In  a  recent  work  we  provided  a  theoretical 
formulation  of  the  problem  in  the  space  and  time  domain 
for  the  reflection-grating  interaction  in  the  nondepleted- 
pump  regime.9  By  using  the  nondepleted-pump  approxi¬ 
mation,  we  reduced  the  two-point  boundary-value  prob¬ 
lem  to  an  initial  value  problem.  In  this  paper  we  present 
a  general  formulation  of  the  problem  in  the  space  and  fre¬ 
quency  domain  based  on  the  standard  statistical  theory 
on  linear  systems.  The  general  formulation  is  also 
implemented  numerically.  Specifically,  we  investigate 
the  signal-intensity  gain  and  the  mutual  coherence  in  the 
contradirectional  wave  mixing  of  two  partially  coherent 
waves.  Contrary  to  conventional  belief,  we  discover  that 
the  effective  interaction  length  (or  grating  length)  can  be 
significantly  longer  than  the  coherence  length  of  the  inci¬ 
dent  waves.  The  results  are  also  compared  with  previous 
theoretical  results  on  simpler  cases  and  with  experimen¬ 
tal  measurements. 


2,  THEORETICAL  MODEL 

Photorefractive  two-wave  mixing  is  a  nonlinear  optical 
process.  Because  of  the  mutual  coherence  of  the  two 
waves,  a  dynamic  holographic  grating  is  formed  in  the 
medium.  Its  position  and  index  profile  are  nearly  sta¬ 
tionary  under  the  condition  of  a  cw  illumination.  Both 
waves  are  scattered  into  each  other  by  the  presence  of 
this  index  grating.  Scattering  of  partially  coherent 
waves  by  a  stationary  grating  can  be  modeled  with  a  sta¬ 
tistical  theory  on  linear  systems.10  An  iterative  proce¬ 
dure  can  subsequently  be  devised  to  obtain  the  final  pho- 
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torefractive  grating  profile  from  an  initially  arbitrary 
grating  profile. 

As  is  shown  in  Fig.  1,  two  counter  propagating  waves 
with  partial  coherence  enter  a  photorefractive  medium  at 
2  =  0  and  z  =  L,  respectively.  The  electric  field  in  the 
photorefractive  medium  can  be  written  as 

E (zf  t)  =  E i(z,  Oexp(-icu0£  +  ik0z) 

+  E2(z ,  Oexp(-io>0^  -  ik0z),  (1) 

where  w0  is  the  center  frequency  of  the  two  partially  co¬ 
herent  waves,  ko  =  nw0/c  is  the  corresponding  wave  vec¬ 
tor,  and  n  is  the  refractive  index  of  the  photorefractive 
medium.  Reflecting  the  partial  coherence,  E^z,  t)  and 
E 2(z,  t)  are  stationary  random  variables.  They  repre¬ 
sent  the  random  fluctuation  of  the  amplitudes  of  the  two 
waves.  For  the  convenience  of  our  later  discussion  we 
will  now  briefly  describe  some  notations  and  definitions 
for  the  second-order  statistical  properties  of  the  two  opti¬ 
cal  waves.  Let  rn(z,  r)  =  (E^z,  ^1)E1*(z,  t2))  and 
r22(z,  r)  =  (E2(z,  ^)E2*(z,  t2))  denote  the  self¬ 
coherence  functions  of  Ex(z,  *x)  and  E^z,  £x),  respec¬ 
tively,  and  r12(z,  r)  =  (Ex(z,  £x)E2*(z,  t2))  be  the 
mutual-coherence  function  between  E^z,  tx)  and 
E2(z,  t2),  where  r  =  tx-  t2  is  the  time  delay  and  () 
means  ensemble  average.  Let  En(z,  A co)  and 
E22(z,  Aw)  denote  the  self-spectral-density  functions  of 
Ex(z,  t)  and  E2(z,  t)f  respectively,  and  E12(zx,  z2,  Aw) 
be  the  cross-spectral-density  function  between  Ex(zx,  t) 
and  E2(z2,  t).  The  spectral-density  functions  and  the 
corresponding  coherence  functions  are  Fourier-transform 
pairs,  i.e., 

Txifz,  r)  =  f  Eu(z,  A w)exp(— iAwr)dAw,  (2) 


r22(^>  r)  =  E22(z,  Aw)exp(-iAwT)dAw, 

J  (3) 


r12(2,  t)  =  Ex2(z,  z,  Aw)exp(-iAwr)dAw, 

J  (4) 


Fig.  1.  Two-wave  mixing  in  photorefractive  crystals  modeled  as 
a  linear  system  with  the  signal-wave  entrance  plane  and  the 
pump-wave  entrance  plane  as  two  input  planes  and  any  arbi¬ 
trary  plane  inbetween  as  the  output  plane. 


where  Aw  =  w  -  w0,  with  w  being  a  general  frequency 
component  of  the  waves.  With  the  above  notations  and 
relations  the  intensity  of  the  two  waves  can  be  expressed 


ikw)dAw, 


7x(z)  ss  ru(z,  0)  =  J  Exx(z,  Aw)dAw,  (5) 

J2(z)  =  r22(z,  0)  =  J  E22(z,  A w)dAw,  (6) 

and  the  mutual  coherence  of  the  two  waves  can  be  ex¬ 
pressed  as 


- 

8.0)  = 

J 


E12(z,  z,  Aw)dAw. 


With  these  equations  the  intensity  of  the  waves  as  well  as 
the  mutual  coherence  can  be  obtained  as  soon  as  the 
spectral-density  functions  are  obtained.  In  what  follows 
we  will  derive  the  spectral-density  functions  by  using  the 
statistical  approach. 

Using  the  above  notations  and  definitions,  we  now  be¬ 
gin  our  discussion  on  the  photorefractive  interaction.  In¬ 
side  the  photorefractive  medium,  a  dynamic  index  grating 
is  generated.  It  can  be  written  as 

y  c  [Q(z,  t) 

Sn  =  -i  - —  T  ,  ,  exp(2ik0z)  +  c.c.  ,  (8) 

2  w0  /0(z) 

where  Q(z,  t)  is  a  measure  of  the  index  grating,  y  is  the 
intensity  coupling  coefficient,  and/0(z)  =  7x(z)  +  /2(z)  is 
the  total  intensity  at  position  z.  For  the  purpose  of  our 
discussion  we  will  call  Ex  the  signal  wave  and  E2  the 
pump  wave.  Thus  for  a  photorefractive  grating  with  a 
positive  %  the  signal  wave  Ex  can  be  amplified.  If  the 
photorefractive  effect  is  based  purely  on  carrier  diffusion 
(e.g.,  BaTi03,)  the  dynamics  of  the  index  grating  is  de¬ 
scribed  by  the  following  equation: 


dQ  (Zy  t) 


+  Q(z,  t)  =  Ex(z,  t)E2*(Zy  t)y 


where  rph  is  the  relaxation-time  constant.  By  virtue  of 
photoexcitations,  photorefractive  processes  are  usually 
slow  at  low  intensities.  It  is  reasonable  to  assume  that 
the  coherence  time  <5oj-1  of  the  two  partially  coherent 
waves  is  much  smaller  than  the  relaxation  time  rph  of  the 
photorefractive  medium,  i.e.,  <5wrph  >  1.  Since  Ex(z,  t) 
and  E2(z,  t)  are  stationary  random  variables,  we  can 
make  the  following  approximation  (see  Appendix  A)8: 

Q(z,  t)  a  (Q(z,  t))  =  r12(z,  0).  (10) 

In  other  words,  two  partially  coherent  waves  with  their 
complex  amplitudes  fluctuating  randomly  with  time  can 
actually  write  a  stationary  grating  in  a  photorefractive 
medium  under  the  appropriate  conditions.  For  simplic¬ 
ity  we  will  denote  Q(z,  t)  with  Q(z)  from  now  on.  Note 
that  Q(z)  is  also  the  mutual  coherence  of  the  two  waves 
at  position  z.  We  note  that  the  approximation  is  also 
valid  when  <5wrph  1. 

Given  arbitrary  functions  of  Q(z)  and  /0(z),  Eq.  (8) 
yields  an  index  grating.  The  propagation  of  a  monochro¬ 
matic  wave  through  such  an  index  grating  can  be  de¬ 
scribed  by  the  coupled  wave  equations.  The  coupled 
monochromatic  waves  can  be  written  as 
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E  =  E i(z,  co)exp (-iojt  +  ikz) 

+  E 2(z,  (o)exp{-io)t  —  ike),  (11) 

where  Ex  and  E2  are  the  amplitudes  of  the  monochro¬ 
matic  components,  o>  is  the  optical  wave  frequency,  and 
k  =  n(i)fc  is  the  optical  wave  vector.  The  coupled  wave 
equations  can  be  written  as 

*Elfe'--  =  27^)  Q(z)®z(z’ 

-  ^E,(2>  «),  (12) 

^2(*’ . —  =  xt4-tQ*(2)Ei(2.  ^)exp(2iAk2) 

<?z  2i0(z) 

+  “E2(z,  o>)#  (13) 

where  Ak  =  k  -  k0  is  the  phase  mismatch  between  the 
optical  waves  and  the  index  grating  and  a  is  the  intensity 
absorption  coefficient.  With  sufficient  boundary  condi¬ 
tions,  Eqs.  (12)  and  (13)  can  be  solved  either  analytically 
in  some  special  cases  or  numerically  in  general.  When 
the  boundary  conditions  are  E x(z  =  0,  w)  =  1  and  E2(z 
=  L,  a))  -  0,  the  solutions  (output)  are  denoted  as 

E^z,  o) )  =  Hn(z,  a) )  and  E 2(z,  co)  =  H12(z,  co).  When 
the  boundary  conditions  are  E j(z  =  0,  co)  =  0  and  E2(z 
=  L,  co)  =  1,  the  solutions  (output)  are  denoted  as 

E^z,  co)  =  H21(z,  co)  and  E2(z,  co)  =  H 22(z,  co).  As  a 
linear  system,  the  general  solutions  are  linear  combina¬ 
tions  of  Hu ,  H12 ,  H21 ,  and  H22 . 

In  each  of  the  iterations  an  arbitrary  stationary  index 
grating  can  be  considered  as  a  linear  system.  Referring 
to  Fig.  1,  we  consider  a  given  stationary  index  grating  in 
the  photorefractive  medium  as  a  linear  system  with  the 
optical  waves  at  the  boundary  planes  z  -  0  and  z  =  L  as 
the  input  and  the  optical  waves  at  an  arbitrary  plane  z  as 
the  output.  According  to  the  theory  on  the  statistical 
properties  of  linear  systems,  the  second-order  statistical 
properties  of  the  optical  waves  at  the  output  plane  can  be 
expressed  in  terms  of  the  second-order  statistical  proper¬ 
ties  of  the  optical  waves  at  the  input  planes  and  the  fre¬ 
quency  response  of  the  linear  system.  To  be  specific,  the 
input  at  the  z  =  0  plane  is  E x(z  =  0,  t),  the  input  at 
z  ~  L  plane  is  E2(z  =  L,  t),  and  the  outputs  at  the  z 
plane  are  Ej(z,  t)  and  E2(z,  t ).  The  frequency  response 
of  this  linear  system  can  be  expressed  in  terms  of  the  so¬ 
lutions  of  Eqs.  (12)  and  (13).  With  the  notations  de¬ 
scribed  in  the  previous  paragraph,  the  frequency  re¬ 
sponses  from  the  input  E^z  =  0,  t)  to  the  outputs 
Ex(z,  t)  and  E2(z,  t)  are 

Hu'(z,  co )  =  Hn(z,  co)exp(iAkz),  (14) 

H12'(z,  co )  =  H12(z,  o>)exp(”iAkz), 

(15) 

respectively.  Similarly,  the  frequency  responses  from 
the  input  E2(z  -  L,  t)  to  the  outputs  Ex(z,  t)  and 
E2(z,  t)  are 

H21'(z,  co )  =  H21(z,  co)exp(iAkL)exp(iAkz),  (16) 


H22'(z,  co)  =  H22(z,  co)exp(iAkL)exp(-iAkz),  (17) 

respectively.  The  additional  phase  terms  in  Eqs.  (14)- 
(17)  account  for  the  difference  Ak  =  k  -  k0  owing  to  a  fi¬ 
nite  A  co  =  co  -  co0.  With  these  spectral-response  func¬ 
tions  we  can  express  the  spectral-density  functions  of  the 
two  outputs  E^z,  t)  and  E2(z,  t)  in  terms  of  the  spectral- 
density  functions  of  the  two  inputs  E x(z  =  0,  t)  and 
E2(z  =  L,  t).  Thus,  according  to  their  definitions  and 
Eqs.  (14)-(17),  we  obtain 

En(z,  Aco)  -  Hn(z,  co)Hu*(z,  co)Eu(z  =  0,  Aco) 

+  H21(z,  co)H21*(z,  co)E22(z  =  L,  Aco) 

+  Hu(z,  co)H21*(z,  co) 

X  E12(z  -  0,  z  =  L,  Aco)exp(-iAkL) 

+  H21(z ,  co)Hn*(z,  co) 

X  [E12(z  =  0,  z  =  L,  Aco) 

X  exp(—  i AkL)]*,  (18) 

E22(z ,  Aco)  =  H12(z,  co)H12*(z,  co)En(z  =  0,  Aco) 

+  H22(z,  co)H22*(z,  co)E22(z  =  L,  Aco) 

+  H12(z,  co)H22*(z,  co) 

X  E12(z  =  0,  z  =  L,  Aco)exp(-iAkL) 

+  H22(z ,  co)H12*(z,  co) 

X  [E12(z  =  0,  z  =  L,  Aco) 

X  exp(-iAkL)]*,  (19) 

E12(z,  z,  Aco)  =  exp(£2Akz){Hu(z,  co)H12*(z,  co) 

X  Eu(z  =  0,  Aco) 

+  H21(z,  co)H22*(z,  co) 

X  E22(z  =  L,  Aco) 

4-  Hn(z,  co)H22*(z,  co) 

X  E12(z  =  0,  z  =  L,  Aco) 

X  exp(-iAkL)  +  H21(z,  co)H12*(z,  co) 

X  [E12(z  =  0,  z  =  Ly  Aco) 

X  exp(-iAkL)]*}.  (20) 

Note  that  the  spectral-density  functions  En(z  =  0,  Aco), 
E22(z  =  L,  Aco),  and  E12(zx  =  0,  z2  =  L,  Aco)  of  the  two 
inputs  E^z  =  0,  t)  and  E2(z  =  L,  t)  are  given  as  the 
boundary  conditions. 

Two  physical  processes  happen  simultaneously  during 
two-wave  mixing  in  a  photorefractive  medium.  First,  the 
two  optical  waves  propagate  through  the  photorefractive 
medium  while  being  scattered  into  each  other  by  the  in¬ 
dex  grating.  Second,  the  scattered  waves  modify  the  in¬ 
dex  grating  through  the  photorefractive  effect  until  a 
steady  state  is  reached.  We  have  provided  above  a  math¬ 
ematical  model  that  describes  these  two  physical  pro¬ 
cesses  separately.  A  steady  state  of  the  two-wave  mixing 
in  the  photorefractive  medium  is  reached  when  the  two 
optical  waves  scattered  by  the  photorefractive  grating  can 
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exactly  sustain  the  same  photorefractive  grating.  A 
steady-state  solution  of  the  two-wave  mixing  in  a  photo- 
refractive  medium  can  thus  be  obtained  by  use  of  the 
mathematical  model  described  above  through  an  iterative 
procedure.  The  procedure  is  outlined  as  follows: 

Step  1:  Give  an  initial  guess  on  the  function 

Q(2)/Z0(z). 

Step  2:  Solve  Eqs.  (12)  and  (13)  for  the  functions 
H fa,  w)  (itj  -  1,2)  with  the  function  Q (z)/J0(z)  pro¬ 
vided  in  the  last  step. 

Step  3:  Obtain  the  spectral-density  functions 
Eii  (z>  A*0)’  E22(z,  Aw),  and  E  12(z,  z,  Aw)  by  use  of  Eqs. 
(18M20). 

Step  4:  Obtain  /^z),  /2(z),  Q(z),  and  Q (z)/IQ(z)  by 
use  of  Eqs.  (5),  (6),  (7),  and  (10). 

Step  5:  Compare  the  new  version  of  the  index  grating 
Q (z)/I0(z)  and  the  previous  version.  If  they  are  close 
within  a  certain  accuracy  requirement,  the  solution  has 
been  obtained.  Otherwise,  the  iteration  continues  by  use 
of  the  new  version  of  the  index  grating. 

To  determine  the  boundary  conditions,  we  assume  that 
both  the  input  optical  waves  are  derived  from  the  same 
laser  source.  If  the  source  laser  wave  has  a  Gaussian 
line  shape  with  a  FWHM  linewidth  of  <$w,  then  the  nor¬ 
malized  spectral-density  function  of  the  source  laser  wave 
can  be  written  as 


4(tt  In  2) 1/2 

Ess(A«>)  =  - - exp 


'  Aw 

2 

— 

l2(ln  2)1/2 

(21) 


Taking  p  as  the  incident-intensity  ratio  Ix{z  -  0)11  %(z 
-  L)  of  the  two  optical  waves  at  their  respective  en¬ 
trance  boundary  planes,  we  obtain  the  boundary  condi¬ 
tions  as 


En(z  -  0,  Aw)  =  /?Ess(Aw),  (22) 

E22(z  =  0,  Aw)  =  Ess(Aw),  (23) 


Ei2(z  1  =  0,  z2  =  Ly  Aw)  =  ^Esa(Aw)exp(“ikoL) 

x  exp(-iw^),  (24) 

where  1 d  is  the  time  delay  between  the  optical  waves 
when  they  reach  their  respective  entrance  planes.  In  de¬ 
riving  the  above  boundary  conditions,  we  have  assumed 
that  the  laser  source  has  a  Gaussian  line  shape.  There  is 
no  loss  of  generality  in  this  assumption.  Similar  results 
can  be  obtained  with  a  different  line  shape. 

This  completes  the  general  formulation  to  model  con- 
tradirectional  two-wave  mixing  in  photorefractive  crys¬ 
tals. 


3.  DISCUSSIONS  AND  SIMULATIONS 

To  clarify  the  complicated  formulation  described  above, 
we  consider  some  simple  cases  first  before  presenting  the 
numerical  and  experimental  results. 

In  the  absence  of  coupling  ( y  =  0 )  in  a  lossless  medium 
(a  =  0)  with  p  =  1  the  above  formulation  describes  the 
interference  of  two  counterpropagating  optical  waves  in  a 


dielectric  medium.  In  this  case  it  takes  only  one  itera¬ 
tion  to  obtain  the  solution,  which  is 

E12(z,  z,  Aw)  =  Ess(Aw)exp(-ikL) 

x  exp(-/wfrf)exp(2i‘Akz).  (25) 

Fourier  transforming  the  above  equation  over  w  and  tak¬ 
ing  into  account  the  extra  phase  term  exp(2ikoz)  resulted 
from  the  definition  of  Ej(z,  t)  and  E2(z,  t)  in  Eq.  (1),  we 
obtain  the  mutual  coherence  of  the  two  waves  as  a  func¬ 
tion  of  position  z , 

I  nL  -  2 nz\ 

r12(z,  r)  =  rj  r  +  td  +  - - -  exp(-iwo^), 

(26) 

where  rs(r)  is  the  normalized  self-coherence  function  of 
the  source  laser  wave  and  n  is  the  index  of  refraction  of 
the  medium.  The  result  has  been  well  established  in  the 
literature.  This  example  can  also  help  us  see  more 
clearly  the  subtle  contribution  of  the  phase  terms  in  the 
above  formulation.  In  our  definition,  td  =  0  ifz  =  L/2  is 
the  plane  of  zero  path  difference. 

When  the  coherence  time  Sio~l  is  much  larger  than  the 
time  delay  involved  in  the  above  formulation,  i.e.,  Sa)td 
<  1  and  ScjnL/c  <  1,  the  normalized  spectral-density 
function  of  the  source  laser  wave  can  be  written  as 
Ess(w)  =  S( w  -  w0).  From  the  above  formulation  we 
can  obtain  the  set  of  equations  governing  the  intensities 
of  the  two  optical  waves, 


d 

hh 

(27) 

8-1 

II 

d 

hh 

(28) 

II 

1-8 

y  +  alx. 

+  12 

These  two  equations  are  exactly  the  same  as  those  ob¬ 
tained  for  two-wave  mixing  of  monochromatic  waves. 
We  note  that  So)td  <  1  is  a  condition  for  partially  coher¬ 
ent  waves  to  be  treated  as  monochromatic  waves  in  con- 
tradirectional  two-wave  mixing  in  a  photorefractive  me¬ 
dium.  We  also  recall  that  the  upper  limit  on  the 
coherence  time,  i.e.,  <5wTph  1,  still  needs  to  be  satisfied 
to  reach  Eqs.  (27)  and  (28)  from  the  above  general  formu¬ 
lation.  This  is  the  main  difference  between  the  results  of 
Eqs.  (27)  and  (28)  in  this  paper  and  those  obtained  di¬ 
rectly  for  monochromatic  waves.  Note  that  we  did  not 
use  the  iterative  procedure  in  obtaining  Eqs.  (27)  and 
(28).  The  complete  boundary  conditions  for  Eqs.  (27)  and 
(28)  are  simply  I\(z  =  0)  and  /2(z  =  L). 

We  now  use  the  above  general  formulation  to  obtain 
the  solution  in  the  nondepleted-pump  regime  (for 
y  =£  0),  which  has  been  obtained  previously  by  a  different 
method.  Using  Eqs.  (12)  and  (13)  and  the  following  three 
definitions, 

Eu(z,  Aw)  =  (Ei(z,  w)E!*(z,  a>))y  (29) 

E22(z,  Aw)  «  (E2(z,  w)E2*(z,  «)>,  (30) 

E12(z,  z,  Aw)  =  <E j(z,  w)E2*(z,  w))exp(2iAkz), 

(31) 

we  can  obtain  a  set  of  equations  governing  the  propaga¬ 
tion  of  the  spectral-density  functions  as 
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<?Eu(z,  Aw) 


—  [QWE12*(2iZ>A,) 

+  Q*(z)Ei2(Z,  Z,  Ao>)] 

—  «En(z,  Aw), 


<?E22(z,  Aw) 


2/0(z) 


[Q(z)E12*(z,  z,  Aw) 


+  Q*(z)E12(z,  z,  Aw)] 
+  <*E22(z,  Aw), 


<7E12(z,  z,  Aw)  _  2iAkEi2(z>  Zj  Aft>)  +  q(2) 

X  [En(z,  Aco)  +  E22(z >  ^Ci))] 

-  aE12(z,  z>  A  <*>).  (34) 

Ei(e,  o>)  and  E2(z,  o>)  in  Eqs.  (29)— (31)  are  related  to  the 
Fourier-transform  coefficients  of  the  two  optical  waves  as 
defined  in  Eq.  (11).  Strictly  speaking,  a  stationary  ran¬ 
dom  process  cannot  be  Fourier  transformed  over  the  time 
variable  t.  However,  we  can  truncate  it  into  a  finite  du¬ 
ration  T  and  Fourier  transform  the  truncated  process. 
After  that,  we  can  first  obtain  a  set  of  equations  similar  to 
Eqs.  (32)-(34)  for  the  truncated  process,  then  let  T  go  to 
infinity  to  obtain  Eqs.  (32)-(34).  This  is  a  standard  pro¬ 
cedure  in  statistical  optics.10  Fourier  transforming  the 
above  relation  over  Aoj,  we  obtain  a  set  of  equations  gov¬ 
erning  the  propagation  of  the  self  and  mutual  coherence 
of  the  two  optical  waves  as 

<?r12(z,  t)  =  2n  dTu(zy  r)  +  y  r12(z,  0) 
dz  c  dr  2  2)  +  I2 

x  [ru(2,  T)  +  r22(*,  r)],  (35) 


dTn(zt  r)  =  T  T12(z9  0) 
dz  2  7i  +  I*> 


Y(Zy  —T) 


,  y  r12*u,  0)  ,  x 

+  2-77T7Tri2(2’r) 


-  «r  u(z,  t), 


<?r22(z,  t)  _  y  r  12(z ,  o) 

<?Z  2  /j  +  1 2 


r12*(z,  -t) 


,  rV(z,0)r  ,  x 

+  2-^77  ri2(z,T) 

+  c*r22(Zy  r).  (37) 

When  there  is  no  absorption  in  the  photorefractive  me¬ 
dium  (i.e.,  a  =  0),  we  can  obtain,  from  Eqs.  (36)  and  (37), 

^[ru(z,  r)  -  T22(z,  r)]  =  0.  (38) 

Further,  letting  r  =  0  in  Eq.  (38),  we  obtain 


yz  Ui(z)  ~  h(z)]  =  0.  (39) 

Equation  (39)  shows  that,  when  there  is  no  absorption  in 
the  photorefractive  medium,  the  intensity  difference  be¬ 


tween  the  signal  wave  and  the  pump  wave  is  a  constant  of 
integration  (conservation  of  power  flow). 

Note  that,  although  Eqs.  (35)-(37)  are  general  and  self- 
consistent,  they  contain  only  the  mutual  coherence  of  the 
two  optical  waves  at  the  same  locations.  They  are  not 
compatible  with  the  kind  of  boundary  conditions  as  given 
by  Eq.  (24),  which  gives  the  mutual  statistical  properties 
of  the  two  waves  at  two  different  points  in  space.  There¬ 
fore  they  are  useful  only  in  the  nondepleted-pump  regime, 
where  we  can  assume  that  the  pump  wave  passes  through 
the  photorefractive  medium  without  changing  its  statisti¬ 
cal  properties.  Under  this  assumption  we  can  obtain 
from  Eqs.  (22)-(24)  the  self-  and  mutual-coherence  func¬ 
tions  of  the  two  waves  at  the  signal- wave  entrance  bound¬ 
ary  plane  (z  =  0)  as 

r12(z  =  0,  t)  =  \lpras(T  +  St)exp(-ia>0St),  (40) 

r„(z  =  0,  r)  ■  pr Jr),  (41) 

r22(z  =  0,  r)  -  r„(r),  (42) 

where  St  =  td  +  nL/c  is  the  time  delay  between  the  two 

optical  waves  at  the  signal-wave  entrance  plane  and 

S(i)T  2 

4(ln  2) m 

is  the  coherence  function  of  the  source  laser  wave.  We 
can  use  Eqs.  (40M42)  as  the  boundary  conditions  and  in¬ 
tegrate  Eqs.  (35)-(37). 

Equations  (35M37)  are  simple  not  only  in  the  sense 
that  they  can  be  implemented  easily  numerically,  but  also 
that  they  can  be  used  to  obtain  approximate  analytical  so¬ 
lutions  to  some  special  cases  within  the  nondepleted  re¬ 
gime.  One  such  case  is  that  the  coherence  length  of  the 
source  laser  wave  is  much  longer  than  the  two-wave¬ 
mixing  interaction  length  and  the  coupling  constant  is 
large.  This  occurs  when  we  use  a  multimode  argon  laser 
and  a  KNb03  :Co  crystal  for  the  two-wave-mixing  experi¬ 
ment.  In  this  case  we  can  neglect  the  term  that  contains 
the  partial  derivative  on  rin  Eq.  (35)  and  reduce  the  set  of 
partial  differential  equations,  i.e.,  Eqs.  (35)-(37),  to  a  set 
of  ordinary  differential  equations  that  can  be  solved  ana¬ 
lytically  under  the  nondepleted-pump  approximation. 
Remember  that  Eq.  (35)  is  derived  directly  from  Eq.  (34). 
We  can  see  from  Eq.  (34)  that  the  approximation  of  ne¬ 
glecting  the  term  that  contains  the  partial  derivative  on  r 
in  Eq.  (35)  implies  that  the  wave-vector  difference  Ak  of 
the  different  frequency  components  of  the  two  partially 
coherent  waves  are  negligible  with  respect  to  the  thick¬ 
ness  of  the  photorefractive  medium  (i.e.,  AkL  ^  1). 
Similar  approximation  has  been  made  previously  by  Sax- 
ena  et  al.  in  the  study  of  multiple-beam  interaction  by 
transmission  gratings  in  the  photorefractive  media.12,13 
Therefore  under  the  nondepleted-pump  approximation 
and  the  approximation  of  AkL  <§  1  the  signal-wave  inten¬ 
sity  gain  and  the  normalized  mutual  coherence  of  the  two 
waves  can  be  obtained  as 

h (z)  __  rn(z,  o)  _  r12(o,  o)r12*(o,  o) 

/i(0)  "  ru(o,o)  rn(o,  o)r22(o,  o) 

X  [exp {yz)  -  ljexp(-az)  +  exp(-az),  (44) 


r8S(r)  =  exp 
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yi2(*>  = 


r  n(z  >  o) 


[r„(*,  o)r22(z,  o)]1/2 


frl2*(o,  o) 
r12(o,  0) 


ru(o,  o)r22(o,  0)  |1/2' 

[1  -  exp(-yz)]  +  -  2— : - exp(-yz)] 

1  12  (u> 


(45) 


In  most  photorefractive  crystals  with  exp[(y-  a)z ]  >  1 
the  signal-intensity  gain  is  affected  primarily  by 
(y  -  a).  We  also  note  that  the  normalized  mutual  co¬ 
herence  of  the  two  waves  is  affected  only  by  the  photore¬ 
fractive  coupling  constant  y. 

As  a  comparison  of  the  general  formulation,  the  simpli¬ 
fied  formulation  and  the  approximate  analytical  solu¬ 
tions,  we  show  in  Figs.  2  and  3  the  signal -intensity  gain 
and  the  normalized  mutual  coherence  of  the  two  waves  at 
the  signal-wave  exit  plane  (z  =  L)  as  a  function  of  the  op¬ 
tical  path  difference  between  the  two  optical  waves  at  the 
signal- wave  entrance  plane  (z  =  0).  The  parameters 
are  y  =  3.0  cm-1,  a  =  0.0  cm"1,  n  —  2.3,  L  =  0.72  cm, 
and  So)  =  2tt  X  1.8  GHz,  and  p  =  10“4.  These  are  typi¬ 
cal  parameters  in  our  experiment  when  we  use  a  multi- 
mode  argon  laser  and  a  KNb03:Co  crystal  to  implement 
the  two-wave-mixing  experiment.  In  this  case  the  coher¬ 
ence  length  of  the  source  laser  wave  is  much  longer  than 
the  thickness  of  the  photorefractive  medium.  The  gen¬ 
eral  formulation  and  the  simplified  formulation  produce 
the  same  results  within  the  numerical  accuracy  in  the 
nondepleted  regime.  The  results  of  these  two  formula¬ 
tions  are  represented  by  the  dashed  curves  in  Figs.  2  and 
3.  We  notice  that  the  approximate  analytical  solution  re¬ 
tains  the  major  characteristics  of  the  exact  solution  and 
provides  a  good  understanding  of  the  interaction.  The 
results  of  the  approximate  analytical  solution  are  repre¬ 
sented  by  the  solid  curves  in  the  figures.  Note  that  the 
curves  obtained  from  the  approximate  analytical  solution 
are  symmetric  about  z0  =  0  owing  to  the  approximation 
of  neglecting  the  partial  derivative  over  r  in  Eq.  (35), 
while  the  curves  obtained  from  the  general  formulation 
are  shifted  to  the  right  side  (see  Figs.  2  and  3).  Here, 
z0  is  defined  as  the  path  difference.  The  difference  be¬ 
tween  the  results  of  the  approximate  analytical  solution 
and  those  of  the  general  formulation  will  decrease  as  the 
ratio  between  the  coherence  length  of  two  waves  and  the 
thickness  of  the  photorefractive  medium  increases. 

Now  let  us  consider  the  case  of  the  depleted  pump. 
Numerical  simulation  with  the  general  formulation  is  the 
only  means  to  analyze  this  case.  We  first  consider  the 
case  in  which  the  coherence  length  of  the  source  laser 
wave  is  finite  but  much  longer  than  the  thickness  of  the 
photorefractive  medium.  Again,  we  use  the  following  set 
of  parameters:  y  —  3.0  cm'1,  a  —  0.0  cm'1,  n  =  2.3, 
L  =  0.72  cm,  and  Su)  =  2rr  X  1.8  GHz.  The  solid 
curves  in  Fig.  4  show  the  signal  intensity  7i(z),  the  pump 
intensity  I2(z),  and  the  normalized  mutual  coherence 
y(z)  =  r12(z,  0)/[ru(z,  0)r22(z,  0)]m  as  functions  of  the 
position  z  in  the  photorefractive  medium.  The  optical 
path  difference  between  the  signal  wave  and  the  pump 
wave  at  the  signal-wave  entrance  boundary  (z  =  0)  is 
chosen  to  be  zero.  The  incidence  intensity  ratio  p  of  the 
two  waves  is  chosen  to  be  one.  We  also  show  in  Fig.  4 
with  dashed  curves  the  signal  intensity  Ix(z)  and  the 


pump  intensity  I2(z)  for  two- wave  mixing  with  monochro¬ 
matic  waves  for  the  purpose  of  comparison.  We  notice 
that  the  results  of  two-wave  mixing  with  partially  coher¬ 
ent  waves  are  very  close  to  those  of  two-wave  mixing  with 


Fig.  2.  Signal-wave  intensity  gain  as  a  function  of  the  optical- 
path  difference  at  the  signal-wave  entrance  plane  in  the 
nondepleted-pump  regime.  The  dashed  curve  is  the  numerical 
solution.  The  solid  curve  is  the  approximate  analytical  solution. 


Fig.  3.  Mutual  coherence  of  the  two  waves  at  the  signal-wave 
exit  plane  as  a  function  of  the  optical  path  difference  at  the 
signal- wave  entrance  plane  in  the  nondepleted-pump  regime. 
The  dashed  curve  is  the  numerical  solution.  The  solid  curve  is 
the  approximate  analytical  solution. 


z(cm) 

Fig.  4.  Signal  intensity,  pump  intensity,  and  the  normalized 
mutual  coherence  as  a  function  of  position  z  in  the  photorefrac¬ 
tive  medium  for  partially  coherent  waves  (solid  curves)  and 
monochromatic  waves  (dashed  curves). 
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monochromatic  waves  when  the  optical  path  difference  of 
the  two  waves  is  small  compared  with  the  coherence 
length  of  the  source  laser  wave.  The  effect  of  partial  co¬ 
herence  of  the  two  waves  in  this  case  is  significant  only 
for  a  large  optical  path  difference  of  the  two  waves.  We 
show  the  signal-intensity  gain  in  Fig.  5(a)  and  the  mutual 
coherence  between  the  signal  wave  and  the  pump  wave  at 
the  pump-wave  entrance  plane  ( z  =  L)  in  Fig.  5(b)  both 
as  functions  of  the  optical  path  difference  between  the  two 
waves  at  the  signal-wave  entrance  plane  (z  =  0)  for  vari¬ 
ous  intensity  ratios  p.  Figure  5(a)  shows  that  the  signal- 
intensity  gain  decreases  as  the  optical  path  difference  of 
the  two  waves  increases  until  there  is  no  coupling  be¬ 
tween  the  two  waves  when  the  optical  path  difference  of 
the  two  waves  exceeds  the  coherent  length  of  the  source 
laser  wave.  The  same  figure  also  shows  that  the  signal- 
intensity  gain  increases  as  the  intensity  ratio  p  decreases 
until  the  signal-intensity  gain  saturates  as  the 
nondepleted-pump  regime  is  reached.  Figure  5(b)  shows 
that  the  normalized  mutual  coherence  of  the  two  waves  at 
the  pump- wave  entrance  plane  (z  =  L)  decreases  as  the 
intensity  ratio  p  of  the  two  waves  increases.  When  p  is 
much  larger  than  one,  the  signal  wave  will  pass  through 
the  photorefractive  medium  with  its  statistical  properties 
almost  unchanged  by  coupling.  In  this  limit  the  normal- 


(a) 


(b) 

Fig.  5.  (a)  Signal-intensity  gain  and  (b)  the  normalized  mutual 
coherence  of  the  two  waves  at  the  pump-wave  entrance  plane 
(z  =  L)  as  functions  of  the  optical  path  difference  at  the  signal- 
wave  entrance  plane  (z  =  0)  for  a  coupling  constant  y 
-  3  cm-1  and  various  intensity  ratios  between  the  signal  wave 
and  the  pump  wave.  Note  that  the  curve  for  a  coupling  constant 
y  -  3  cm-1  and  ft  >  1  is  the  same  as  that  for  a  coupling  constant 
7=0  cm-1  and  an  arbitrary  p. 


ized  mutual  coherence  of  the  two  waves  at  the  pump-wave 
entrance  plane  (z  =  L)  with  coupling  will  be  the  same  as 
that  for  an  arbitrary  p  but  without  coupling.  We  note 
that  the  normalized  mutual  coherence  of  the  two  waves  at 
the  pump-wave  entrance  plane  (z  =  L)  is  enhanced  by 
coupling.  Figure  5(b)  also  shows  that  the  normalized 
mutual  coherence  decreases  quickly  as  the  optical  path 
difference  gets  close  to  and  larger  than  the  coherence 
length  of  the  source  laser  wave. 

Second,  we  consider  the  case  that  the  coherence  length 
of  the  source  laser  wave  is  shorter  than  the  thickness  of 
the  photorefractive  medium.  The  following  parameters 
are  chosen  in  the  simulation:  a  =  0.0  cm-1,  n  =  2.3, 
L  =  2.0  cm,  S<o  =  2  ir  X  18  GHz,  and  td  =  - nL/c .  td 
=  -nLIc  implies  that  the  optical  path  difference  be¬ 
tween  the  signal  wave  and  the  pump  wave  is  zero  at  the 
signal-wave  entrance  plane  ( z  =  0).  Since  the  spectral 
line  shape  of  the  source  laser  wave  is  assumed  to  be 
Gaussian,  the  coherence  length  of  the  source  laser  wave 
inside  the  photorefractive  medium  is  Lc  =  2tt 
x  0.664c/(rt££u)  =  0.48  cm.  In  Fig.  6  we  show  the  sig¬ 
nal  intensity  Ii(z)  and  the  grating  profile  Q (z)/J0(z)  as 
functions  of  the  position  z  inside  the  photorefractive  me¬ 
dium  for  an  incident  intensity  ratio  p  =  1  and  a  coupling 
constant  y  =  20  cm-1.  We  note  that  the  length  of  the 
photorefractive  grating  is  limited  by  the  partial  coherence 
of  the  two  interacting  waves  and  that  the  pump  depletion 
is  moderate  even  for  a  very  large  coupling  constant.  In 
Figs.  7(a)  and  7(b)  we  show  the  signal  intensity  /a(z)  and 
the  grating  profile  Q (z)//0(z)  as  functions  of  the  position 
z  inside  the  photorefractive  medium  for  an  incident  inten¬ 
sity  ratio  p  =  10"4  and  coupling  constants  y  -  10  cm-1 
and  y  -  20  cm"1,  respectively.  We  note  that  the  length 
of  the  photorefractive  grating  is  increased  but  still  prima¬ 
rily  limited  by  the  partial  coherence  of  the  two  interacting 
waves  for  a  small  incident  intensity  ratio  and  a  small  cou¬ 
pling  constant.  Figure  7(b)  shows  that  the  length  of  the 
photorefractive  grating  is  no  longer  limited  by  the  partial 
coherence  of  the  two  interacting  waves  for  a  small  inci¬ 
dent  intensity  ratio  and  a  large  coupling  constant.  The 
length  of  the  photorefractive  grating  in  this  case  is  lim¬ 
ited  by  the  length  of  the  photorefractive  medium,  which  is 
much  longer  than  the  coherence  length  of  the  interacting 
waves.  The  photorefractive  grating  is  a  temporally  sta¬ 
tionary  index  grating.  When  the  incident  intensity  ratio 
is  small,  the  amplified  signal  wave  at  its  exit  plane  (z 
=  L)  is  primarily  the  incident  pump  wave  reflected  by 
the  photorefractive  grating,  and  the  depleted  pump  wave 
at  its  exit  plane  (z  =  0)  is  primarily  the  incident  pump 
wave  transmitted  through  the  photorefractive  grating. 
When  the  length  of  the  photorefractive  grating  is  compa¬ 
rable  to  or  longer  than  the  coherence  length  of  the  source 
laser  wave  inside  the  photorefractive  medium,  the  output 
waves  Ej(L)  and  E2(0)  may  have  spectra  different  from 
those  of  the  input  waves  E^O)  and  E2(L)  owing  to  the 
presence  of  the  photorefractive  grating,  which  acts  as  a 
spectral  filter.  In  Figs.  8(a)  and  8(b)  we  show  the  nor¬ 
malized  spectra  of  the  amplified  signal  wave  and  the  de¬ 
pleted  pump  wave  at  their  respective  exit  planes,  using 
the  same  sets  of  parameters  as  those  used  to  obtain  Figs. 
7(a)  and  7(b),  respectively.  We  note  that  the  amplified 
signal  wave  has  a  bandwidth  narrower  than  that  of  the 
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Fig.  6.  Signal  intensity  Ix(z)  (solid  curve)  and  the  grating  pro¬ 
file  Q (z)ll0(z)  (dashed  curve)  as  functions  of  the  position  z  inside 
the  photorefractive  medium  for  an  incident  intensity  ratio  (3 
=  1  and  a  coupling  constant  y  =  20  cm-1. 
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pump  wave  has  a  different  spectrum  from  the  incident 
pump  wave.  According  to  the  spectra  shown  in  Fig.  8, 
the  case  considered  in  Fig.  7(a)  and  Fig.  8(a)  is  in  the 
nondepleted-pump  regime,  and  the  case  considered  in 
Fig.  7(b)  and  Fig.  8(b)  is  in  the  depleted-pump  regime,  al¬ 
though  both  cases  seem  to  be  in  the  nondepleted-pump  re¬ 
gime  when  we  look  only  at  the  signal  intensity.  In  Figs. 
9(a)  and  9(b)  we  show  with  solid  curves  the  normalized 
mutual  coherence  of  the  two  waves  as  a  function  of  posi¬ 
tion  z ,  again  using  the  same  sets  of  parameters  as  those 
used  to  obtain  Figs.  7(a)  and  7(b),  respectively.  For  com¬ 
parison  we  show  in  these  two  figures  with  dashed  curves 
the  normalized  mutual  coherence  of  the  two  waves  as  a 
function  of  position  z  without  coupling.  In  Fig.  9(a)  the 
normalized  mutual  coherence  of  the  two  waves  is  the 
same  with  coupling  as  that  without  coupling  in  the  region 
z  <  0  and  is  increased  in  the  region  z  >  L  owing  to  cou¬ 
pling;  in  Fig.  9(b)  the  normalized  mutual  coherence  of  the 
two  waves  is  decreased  in  the  region  close  to  plane 
2  =  0.  The  increase  of  mutual  coherence  in  both  Figs. 
9(a)  and  9(b)  can  be  attributed  to  the  reflection  of  the 
strong  pump  wave  in  the  direction  of  the  weak  signal 
wave  by  the  stationary  photorefractive  index  grating. 
The  decrease  of  mutual  coherence  in  Fig.  9(b)  can  be  at¬ 
tributed  to  the  spectral-filtering  effect  of  the  photorefrac¬ 
tive  index  grating  on  the  pump  wave.  There  is  no  de¬ 
crease  of  mutual  coherence  in  the  region  close  to  plane  z 
=  0  in  Fig.  9(a)  because  pump  depletion  and  therefore 
the  spectral-filtering  effect  on  the  pump  wave  are  negli¬ 
gible  in  this  case. 

Third,  when  a  laser  beam  enters  a  photorefractive  crys¬ 
tal,  scattering  occurs  because  of  surface  pits,  imperfec- 


(b) 


Fig.  7.  Signal  intensity  Ix{z)  (solid  curve)  and  the  grating  pro¬ 
file  Q  (z)//0(z)  (dashed  curve)  as  functions  of  the  position  z  inside 
the  photorefractive  medium  for  an  incident  intensity  ratio  /3 
=  1CT4  and  coupling  constants  (a)  y  =  10  cm"1  and  (b)  y 
=  20  cm"1. 

incident  pump  wave.  Furthermore,  the  bandwidth  of  the 
amplified  signal  wave  decreases  as  the  length  of  the  pho¬ 
torefractive  grating  increases.  In  Fig.  8(a),  since  the  cou¬ 
pling  constant  is  small,  the  pump-wave  depletion  is  small 
( y  =  10  cm-1)  and  the  spectrum  of  the  transmitted  pump 
wave  is  almost  the  same  as  that  of  the  incident  pump 
wave.  In  Fig.  8(b),  since  the  coupling  constant  is  large 
(y  =  20  cm"1),  the  pump  depletion  is  significant.  The 
central  part  of  the  incident  pump-wave  spectrum  is  de¬ 
pleted  most  significantly,  and  therefore  the  transmitted 
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Fig.  8.  Normalized  spectra  of  the  amplified  signal  wave  (solid 
curve)  and  the  depleted  pump  wave  (dashed  curve)  at  their  re¬ 
spective  exit  planes  for  an  incident  intensity  ratio  an(l 

coupling  constants  (a)  y  =  10  cm" 1  and  (b)  y  =  20  cm  . 
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Fig.  9.  Normalized  mutual  coherence  of  the  two  waves  as  func¬ 
tions  of  position  z  for  an  incident  intensity  ratio  (3  =  10  "4  and  for 
coupling  constants  (a)  y  =  10  cm-1  (solid  curve)  and  y 
=  0  crrT 1  (dashed  curve)  and  (b)  y  =  20  cm'1  (solid  curve)  and 
y  =  0  cm-1  (dashed  curve). 


Fig.  10.  Phase-conjugation  reflectivity  as  a  function  of  the  co¬ 
herence  length  of  the  incident  beam  for  various  value  of  the  cou¬ 
pling  constant  y. 

tion,  and  defects  in  the  crystal.  The  scattered  light  over¬ 
laps  with  the  incident  beam  and  they  undergo  two-wave 
mixing.  Under  the  appropriate  condition  the  scattered 
light  can  be  amplified,  leading  to  phenomena  such  as  fan¬ 
ning  and  stimulated  backscattering.  In  a  manner  very 
similar  to  stimulated  Brillouin  scattering  the  stimulated 
backward  scattering  in  photorefractive  media  is  a  pos¬ 
sible  mechanism  for  self-pumped  phase  conjugation.14,15 
The  general  formulation  developed  in  this  paper  can  be 
employed  to  investigate  the  effect  of  coherence  on  self- 
pumped  phase  conjugation  by  2k  gratings.  We  show  in 
Fig.  10  the  phase-conjugation  reflectivity  as  a  function  of 


the  coherence  length  of  the  incident  beams  for  various 
value  of  the  coupling  constant  y.  The  parameters  in  this 
simulation  are  n  -  2.3,  L  =  0.72  cm,  a  =  0.0  cm"1,  and 
f3  =  1.0  x  10“6.  The  time  delay  of  the  two  waves  is  as¬ 
sumed  to  be  zero  at  the  signal-wave  entrance  plane 
( z  —  0).  The  coherence  length  of  the  incident  beams  is 
related  to  the  bandwidth  as  Lc  =  2tt  X  0.664 c/Sw.  We 
note  that  the  phase-conjugation  reflectivity  increases  as 
the  coherence  length  of  the  incident  beam  increases  and 
reaches  to  a  constant  when  the  coherence  length  of  the  in¬ 
cident  beam  is  much  longer  than  the  thickness  of  the  pho¬ 
torefractive  medium.  We  also  note  that  the  phase- 
conjugation  reflectivity  increases  as  the  coupling  constant 
increases. 


4.  EXPERIMENTS 

The  above  theory  is  validated  experimentally.  The  ex¬ 
perimental  setup  is  shown  in  Fig.  11.  We  utilized  a  45°- 
cut  KNb03:Co  crystal  (the  c  axis  is  in  the  horizontal 
plane  leaning  toward  the  z  =  0  face  of  the  crystal,  and 
the  b  axis  is  in  the  vertical  direction).  The  measured  pa¬ 
rameters  of  the  crystal  are  y  ~  3.3  cm-1,  a  =  0.5  cm"1, 
n  =  2.3,  and  L  =  0.72  cm.  A  multimode  argon  laser  op¬ 
erating  at  514  nm  is  used  as  the  laser  source  with  a  mea¬ 
sured  FWHM  bandwidth  of  1.83  GHz.  The  extraordi¬ 
nary  polarization  of  the  laser  wave  is  used  in  the 
experiment.  As  illustrated  in  Fig.  11,  the  signal  wave 
and  the  pump  wave,  obtained  by  splitting  the  argon  laser 
wave,  propagate  contradirectionally  into  the  KNb03:Co 
crystal.  The  incident  intensity  ratio  of  the  two  waves  is 
p  -  0.00151.  The  power  of  the  pump  wave  is  main¬ 
tained  at  —50  mW.  The  optical  path  difference  of  the 
two  waves  at  the  signal-wave  incident  plane  z  =  0  is  de¬ 
noted  as  AL  =  L2  -  Lx .  To  monitor  the  mutual  coher¬ 
ence  between  the  signal  wave  and  the  pump  wave  at  the 
output  plane  z  =  L,  we  employed  another  reference  wave 
(E2ref)  that  was  split  from  the  pump  wave  E2 .  The  opti¬ 
cal  path  difference  of  Ex  and  E2ref  waves  was  adjusted  to 
be  the  same  as  that  of  Ex  and  E2  waves  at  the  output 
plane  z  =  L.  By  a  simple  homodyne  technique,  the  in- 


Fig.  11.  Experimental  setup.  The  distances  Lx  and  L2  are  the 
optical  path  length  of  the  signal  wave  and  the  pump  wave  from 
laser  source  to  the  signal-wave  incident  plane  2  =  0,  respec¬ 
tively.  L2ref  is  the  optical  path  length  of  reference  wave  Z£2re f 
from  the  laser  source  to  the  signal  output  plane  z  =  L. 
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(a)  (b) 

Fig.  12.  Interference  pattern  of  the  signal  wave  and  the  refer¬ 
ence  wave  at  the  output  plane  PI  (a)  without  photorefractive  cou¬ 
pling  and  (b)  with  photorefractive  coupling.  Note  the  increase  of 
fringe  visibility  that  is  due  to  the  coupling. 


are  also  provided  for  various  coupling  constants  and  vari¬ 
ous  incident  intensity  ratios  between  the  signal  wave  and 
the  pump  wave.  We  found  that  the  mutual  coherence  be¬ 
tween  the  signal  wave  and  the  pump  wave  can  be  both  in¬ 
creased  and  decreased  owing  to  coupling.  We  also  found 
that  both  the  strength  and  the  length  of  the 
photorefractive-index  grating  increases  as  the  coupling 
constant  increases.  A  photorefractive  index  grating 
much  longer  than  the  coherence  length  of  the  incident 
waves  can  be  formed  when  the  coupling  constant  is  large 
and  the  incident  intensity  ratio  is  small.  Owing  to  the 
spectral-filtering  effect  of  the  photorefractive-index  grat¬ 
ing,  the  spectra  of  the  two  interacting  waves  will  be  al¬ 
tered  as  they  pass  through  the  photorefractive  medium. 
The  general  formulation  is  further  used  to  model  the  ef¬ 
fect  of  partial  coherence  on  self-pumped  phase  conjuga¬ 
tion  by  a  2k  grating.  The  theoretical  predictions  are  in 
excellent  agreement  with  experimental  measurements. 


optical  path  difference  (cm) 

Fig.  13.  Signal-intensity  gain  as  a  function  of  the  optical  path 
difference  of  the  two  waves  at  the  signal-wave  entrance  plane 
( z  =  0).  The  dots  are  experimental  data,  and  the  solid  curve  is 
the  theoretical  data. 

terference  fringes  generated  by  Ex  and  E2ref  waves  were 
observed  by  a  CCD  camera  at  the  output  plane  PI.  The 
normalized  mutual  coherence  y12(L)  can  be  estimated  as 
(7max  -  where  (7max  -  7^)  is  the  ampli¬ 

tude  of  the  fringes.  In  our  experiment  we  monitored  the 
interference  pattern  with  and  without  pump  beam  E2. 
Figure  12  shows  the  interference  patterns  with  a  normal¬ 
ized  mutual  coherence  r12(0,  0)  0.43  at  z  =  0  (AL 

=  4  cm).  The  measured  normalized  mutual  coherence 
increases  from  0.19  to  0.7  at  z  =  L.  The  intensity  gain 
of  the  signal  wave  was  also  measured.  Figure  13  shows 
the  measurement  of  the  intensity  gain  (dots)  of  the  signal 
wave  at  the  z  =  L  plane  as  a  function  of  the  optical  path 
difference  AL.  Along  with  the  data  is  the  theoretical 
curve  for  the  same  parameters.  An  excellent  agreement 
between  theory  and  experiment  was  achieved. 

5.  CONCLUSIONS 

In  conclusion,  we  have  investigated  theoretically  contra- 
directional  two-wave  mixing  with  partially  coherent 
waves  in  photorefractive  crystals.  A  general  formulation 
based  on  the  theory  of  statistical  properties  of  linear  sys¬ 
tems  is  provided.16  Previous  results  on  several  simpli¬ 
fied  cases  are  rederived  as  special  cases  of  the  general  for¬ 
mulation  so  that  we  can  get  more  insight  into  the  general 
formulation  as  well  as  into  the  simplified  cases.  Results 
of  numerical  implementation  of  the  general  formulation 


APPENDIX  A:  ERGODICITY  PROPERTY 

In  this  appendix  we  provide  the  derivation  of  Eq.  (10). 
Equation  (9)  gives  the  dynamics  of  the  photorefractive  in¬ 
dex  grating.  It  can  be  rewritten  into  an  integration  form 
as 

Q(2>  t)  =  L  J  EjCz,  OE2*(2,  r)exp(~^~)dr- 

Note  that  exp[(f'  -  <)/Tph]dt'  =  1.  According  to 

Eq.  (Al),  the  grating  amplitude  Q(z,  t)  is  approximately 
the  average  value  of  Ei(z,  t') E2*(z,  t')  over  a  time  pe¬ 
riod  rph . 

Since  the  optical  wave  amplitudes  Ex(z,  t)  and 
E2(z,  t)  are  stationary  random  processes,  the  ensemble 
average  (Ei(z,  £)E2*(z,  £))  is  independent  of  the  time 
variable  t .  Taking  the  ensemble  average  of  Eq.  (Al),  we 
can  obtain 

(Q(z,  t))  -  <Ex(z,  *)E2*(z,  t))  -  r12(z,  0).  (A2) 

Equation  (A2)  shows  that  the  ensemble  average  of  the 
grating  amplitude  Q(z,  t)  is  equal  to  the  mutual  coher¬ 
ence  of  the  two  waves. 

In  general,  the  grating  amplitude  Q(z,  t)  is  also  a  ran¬ 
dom  variable.  It  fluctuates  around  its  ensemble  average. 
The  mean  square  value  of  this  random  fluctuation,  also 
called  the  variance  of  the  random  variable  Q(z,  t ),  can  be 
written  as 

<|Q(Z,  t)  -  <0(2,  *)>l2>  =  (10(2,  <)l2>  -  1(0(2,  <)>la. 

(A3) 

By  using  Eq.  (Al),  we  can  write  the  first  term  on  the  right 
side  of  Eq.  (A3)  explicitly  as 

(|Q(z,0l2)=  A  f  f  <FU, 

Tph  J  —CO  J  —CO 

x  expjf— ^)exp(^^)di'd«'', 

(A4) 
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where  F(z,  t)  =  E t(z,  t)E2*(zt  t)  is  a  shorthand  nota¬ 
tion.  Since  the  optical  wave  amplitudes  Ex(z,  t)  and 
E 2(z,  t)  are  stationary  random  processes,  F (z,  t)  is  also  a 
stationary  random  process.  Therefore  the  ensemble  av¬ 
erage  (F (z,  t')F*(z,  t"))  is  a  function  of  only  two  vari¬ 
ables  2  and  tr  ~  t”.  We  define  R(z,  V  -  t”) 
=  (F(z,  7')F*(z,  tM))  as  a  shorthand  notation.  With 
this  relation  we  can  change  the  integration  arguments  in 
Eq.  (A4)  from  t'  and  t"  to  t\  —  t'  +  tn  —  2 1  and  t2  =  t* 
-  t"  and  rewrite  Eq.  (A4)  as 

(|Q(2,  ol2)  =  r~2  f  d<i  exP  — )  f  d<2R(z,  t2). 

"Tph  J \  ^ph /  J 

(A5) 

Note  that  l/(2Tph2)/°„d<1  exp(i1/Tph)/‘<1d<2  =  1.  Sub¬ 
stituting  Eqs.  (A2)  and  (A5)  into  Eq.  (A3),  we  can  rewrite 
the  variance  of  the  grating  amplitude  Q (2,  t )  as 

<IQ(z,  t)  -  (Q (z,  t))|2> 

/,>«*•« 

-  |r12(z,  0)|2].  (A6) 

Remember  that  F(z,  t)  is  a  stationary  random  process. 
Let  A  7  be  the  minimum  time  delay  that  is  necessary  for 
F(z,  t)  and  F(z,  t  ±  A7)  to  be  uncorrelated.  For  \t2\ 
^  A7,  R(z,  t2)  -  |r12(z,  0)|2  is  equal  to  zero.  For  jf2| 
^  A 7,  R (2,  t2)  -  |r12(2,  0)|2  is  a  function  of  t2.  Since 
both  R(z,  72)  and  |r12(z,  0)|2  are  of  the  order  of 
71(2)72(2),  the  upper  bound  of  R(z,  t2)  -  |r12(2,  0)|2  can 
be  written  as  m71(2)72(2),  where  m  is  a  constant  factor  of 
the  order  of  unity.  With  these  estimations  we  can  obtain 
from  Eq.  (A6) 

<|Q (2,  t)  -  (Q (2,  0)|2)  <  mli(z)l2{z)  — .  (A7) 

'  ph 

According  to  relation  (A7),  the  fluctuation  of  the  grating 
amplitude  Q(z,  t)  decreases  as  the  ratio  A £/rph  decreases. 
When  A  t/rvh<  1,  we  can  neglect  the  fluctuation  of  the 
grating  amplitude  Q (2,  t)  and  obtain  Eq.  (10),  i.e., 

Q (2,  0  -  (Q (2,  O)  =  r12(z,0).  (A8) 

Usually  A t  is  of  the  order  of  the  coherence  time  of 

the  incident  waves.  In  this  case  the  condition  for  Eq. 
(A8)  can  also  be  written  as  <5wrph  1. 

In  this  appendix  we  have  shown  that  the  time  average 
of  the  random  variable  E^z,  OE2*(z,  t)  is  equal  to  the 
ensemble  average  of  the  random  variable 


Ex(z,  OE2*(z,  t).  This  type  of  property  is  referred  to  as 
ergodicity  in  statistics. 
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